Naval  Research  Laboratory 

Washington,  DC  20375-5320 


AD-A266  105 


NRL/MR/6650--93-7357 


Theory  and  Modeling  of 
Thermal  Reaction  Propagation  in 
Beam-Initiated  Explosives 


J.  B.  Aviles,  Jr. 
Consultant,  SFA,  Inc . 


A.  Stolovy 

Directed  Energy  Effects  Branch 

Condensed  Matter  and  Radiation  Sciences  Division 


June  10, 1993 


93 


■£> 


93-14022 

Ml  till!  *»<  .  J 


I  N: 


M) 


Approved  for  public  release,  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  Information  ie  estimated  to  avenge  1  hour  par  reaponte.  Including  the  time  for  reviewing  Instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  Information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  Information,  Including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1216  Jefferson 

Davis  Highway,  Suite  1 204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-01 88).  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  {Leave  Blank) 

i.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

June  10.  1993 

Final  Report 

4.  TITLE  ANO  SUBTITLE 

5.  FUNDING  NUMBERS 

Theory  and  Modeling  of  Thermal  Reaction  Propagation  in  Beam-Initiated  Explosives 

PE  -  61153N-12 

TA  -  RRO 12-02-44 

6.  AUTHOR(S) 

WU  -  66-2821-0-3 

J.B.  Aviles  Jr.*  and  A.  Stolovy 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSIES) 

Naval  Research  Laboratory 

Washington.  DC  20375-5320 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRUMR/6650-93-7357 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

1 1 ,  SUPPLEMENTARY  NOTES 

"Consultant,  SFA,  Inc. 

12*.  DISTRIBUTION/A VAILABIUTY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 

1 1 3.  ABSTRACT  Maximum  200  word*) 

Previous  work  on  the  thermal  initiation  of  energetic  materials  by  particle  beams  concentrated  on  uniformly  irradiated  confined 
samples.  The  present  report  extends  the  work  to  situations  where  only  a  portion  of  the  material  has  been  irradiated.  In  some 
materials  (e.g.,  TATB)  only  the  irradiated  region  undergoes  thermal  Initiation;  in  other  materials  a  thermal  initiation  front  has 
been  observed  (by  thermocouple  response)  to  propagate  away  from  the  irradiated  region.  The  frontal  speeds  are  of  the  order  of 
centimeters  per  second  to  tens  of  centimeters  per  second.  These  slow  speeds  eliminate  the  possibility  of  pressure  propagation  or  a 
deflagration  to  detonation  transition  (DDT), 

Modeling  is  based  on  the  interference  that  at  the  thermal  initiation  temperature  the  material  develops  a  greatly  enhance 
thermal  and  gaseous  diffusivity.  This  allows  the  exothermic  energy  reteased  to  propagate  into  the  cold  region  and  forestalls  any 
explosive  temperature  rise.  At  a  well  defined  positive  R*  (the  explosion  radius)  thermal  initiation  propagation  ceases  and  an 
explosive  temperature  rise  takes  place. 

14.  SUBJECT  TERMS 

16  NUMBER  OF  PAGES 

63 

10.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

18.  SECURITY  CLASSIFICATION 

OF  THIS  RAGE 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  7640-01*280-6800 

SuMvdfwm  ill  tft«v.  2-8)1 

A*»tf*W6V  AKSl  5IJ2MU 


m.icj 


CONTENTS 


INTRODUCTION  AND  SUMMARY .  1 

1.  EXPERIMENTAL  OBSERVATION  .  4 

2.  CONCEPTUAL  PHYSICAL  DESCRIPTION  .  6 

3.  EXOTHERMAL  REACTION  MODELING .  10 

4.  ADIABATIC  THERMAL  INITIATION  BY  PARTICLE  BEAMS .  13 

5.  SIMPLE  MODEL  FOR  INITIATION  PROPAGATION  .  16 

6.  GENERAL  TRANSPORT  THEORY  AND  MODELING  . . .  20 

7.  SOLUTION  OP  THE  MODEL  EQUATION  .  38 

8.  MODEL  PARAMETERS  FROM  EXPERIMENTS  .  46 

9.  CONCLUSION . 50 

APPENDIX  A.l.  Early  NRL  Propagation  Experiment .  52 

APPENDIX  A.2.  Reaction  Model  Equivalencies  .  57 

REFERENCES  .  62 

ACKNOWLEDGMENT  .  63 


ioooasion  For 

HITS  GRA&I 

ar 

DUG  TAB 

G 

Uuuunocaood 

O 

JuatlftC;aloa. — 

Ww  _ -J 

Distribution/ _ 

»  •  *  1 

Availability  Cofloa 
Avail 

Blot  Sp««U'l 


THEORY  AND  MODELING  OF  THERMAL  REACTION  PROPAGATION 

IN  BEAM-INITIATED  EXPLOSIVES 

Introduction  and  Summary 

Over  the  past  years  the  thermal  initiation  properties  of  energetic  materials 
have  been  studied  by  Stolovy,  Namenson,  Kidd,  Jones,  and  Aviles  [Ref.  1].  These 
studies  have  overwhelmingly  concentrated  on  the  behavior  of  small  confined 
samples  uniformly  irradiated  by  particles  beams.  Since  thermal  losses  are  greatly 
reduced,  the  thermal  behavior  is  essentially  adiabatic.  The  resultant  thermal 
initiation  threshold  energies  and  temperatures  are  important  intrinsic  properties 
of  the  material.  In  particular  there  is  direct  applicability  to  charged  and  neutral 
particle  beam  weapons  lethality  in  scenarios  where  energetic  materials  are 
irradiated  in  depth. 

This  report  is  concerned  with  effects  that  occur  after  initiation  in  a 
situation  where  only  a  portion  of  the  material  has  been  irradiated.  Given  the 
rapidity  with  which  the  temperature  rises  after  initiation  in  the  adiabatic  case 
and  given  the  small  thermal  conductivity  of  energetic  materials,  it  would  be 
natural  to  predict  that  thermal  initiation  will  occur  only  in  the  beam  irradiated 
region  that  has  reached  the  thermal  initiation  threshold.  However,  for  certain 
energetic  material,  Stolovy,  et  al  have  observed  thermal  initiation  fronts 
propagating  away  from  the  beam  region.  The  frontal  speeds  of  various  materials 
are  of  the  order  of  centimeters  per  second  to  tens  of  centimeters  per  second.  The 
relatively  low  speeds  eliminate  the  possibility  of  pressure  waves  and  the  light 
confinement  is  not  conducive  to  a  deflagration  to  detonation  transition  (DDT). 
Thus  some  kind  of  thermal  reaction  propagation  mechanism  is  indicated.  The 
measured  thermal  diffusivities  at  subinitiation  temperatures  are  too  low  by 
orders  of  magnitude  to  account  for  these  speeds.  It  is  necessary  to  postulate  that 
a  radical  transition  in  transport  properties  occurs  at  the  thermal  initiation 
temperature.  Model  equations  have  been  developed  for  this  situation  and  have 
been  solved  analytically  (in  an  appropriate  approximation)  for  the  cylindrically 
symmetric  case  which  obtains  in  the  electron  beam  experiments.  In  particular, 
the  time  for  the  initiation  front  to  propagate  from  the  edge  of  the  particle  beam  of 
radius  a  to  the  radial  position  R  is  given  by 

At  -  tiVl-(R/R’)2+ to  (wHR/R  )2ff  (0.1) 

where 


R'=f2Dft, 


(0.2) 


The  parameters  tj  is  the  time  for  the  temperature  in  the  cold  high  explosive  (HE) 
to  rise  to  *he  initiation  temperature  as  a  result  of  propagation.  The  parameter  t*  is 
the  adiabatic  time  to  explosion,  i.e.,  the  time  after  initiation  for  the  HE  to  reach  a 
very  high  temperature  when  there  is  no  outflow  of  heat.  The  parameter  Dr  is  an 
effective  thermal  diffusivity  which  pertains  only  to  region  of  the  HE  which  is 
above  the  thermal  initiation  temperature. 
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The  most  important  feature  of  Eq.  (1)  is  that  for  values  of  R  >  R*  the 
expression  takes  on  complex  values.  The  theoretical  analysis  shows  that  when 
the  front  reaches  R*  the  initiated  region  R  <  R*  undergoes  a  sudden  rapid  rise  in 
temperature,  i.e.,  an  explosion  develops.  Hence,  R*  will  be  called  the  explosion 
radius  for  the  material.  If  the  beam  radius  exceeds  the  explosion  radius  no 
propagation  occurs;  the  irradiates  region  then  undergoes  an  adiabatic  thermal 
explosion. 

The  experimental  data  for  the  validation  of  Eqs.  (0.1)  and  (0.2)  are  sparse. 
The  initial  experiments,  which  discovered  the  phenomenon  of  thermal  initiation 
propagation  in  HE,  were  of  an  exploratory  nature.  The  results,  not  previously 
published,  are  described  in  Appendix  A.l.  Unfortunately,  the  set  of  experiments 
[Ref.  2]  which  were  designed  to  investigate  the  phenomenon  in  greater  detail  did 
not  have  the  benefit  of  the  present  analysis.  The  beam  radii  were  generally 
chosen  larger  than  in  the  initial  experiments  and  the  thermocouples  used  were 
not  as  responsive.  As  a  result  propagation  was  obscured.  These  later  experiments 
are  compatible  with  the  theory  presented  here  but  do  not  contribute  to  its 
quantitative  validation.  The  initial  experiments  on  FBX-9404,  although  sparse  in 
data,  are  sufficient  to  give  order  of  magnitude  values  for  the  parameters  ti  and 
R*.  An  estimate  of  t,  from  thermal  decomposition  reaction  parameters  then 
allows  a  determination  of  Dr. 

Subsequent  experiments,  designed  to  validate  the  present  theory,  have  not 
been  carried  out  owing  to  funding  limitations.  Rather  than  allow  the 
phenomenon  of  thermal  initiation  propagation  to  remain  obscure  and 
unexplored,  a  fairly  ambitions  conceptual  and  theoretical  analysis  was 
undertaken.  The  genera)  theory  and  the  results  for  particle  beam  irradiation,  i.e., 
for  cylindrical  symmetry,  are  present  here,  laser  irradiation  of  a  surface  (planar 
symmetry)  also  produces  thermal  initiation  propagation  and  has  been  analyzed 
theoretically;  however  the  results  are  not  present  here  since  there  is  insufficient 
data  to  make  a  quantitative  analysis.  These  experiments  are  also  described  in 
Appendix  A.!;  they  clearly  indicate  a  propagating  thermal  initiation  front 

In  Section  1  the  experimental  observations  and  certain  reasonable 
conclusions  are  given.  Details  are  found  in  Appendix  A.1  and  Reference  2.  In 
Section  2  a  conceptual  physical  description  is  presented;  it  is  based  on  the 
experimental  observations,  the  model  which  is  developed  in  later  sections,  and 
reasonable  speculation.  Empirical  reaction  modeling,  based  on  thermal 
decomposition  experiments,  is  presented  in  Section  3  and  is  applied  in  Section  4 
to  thermal  initiation  of  materials  uniformly  irradiated  by  particle  beams. 
Appendix  A.2  shows  that,  for  the  limited  temperature  range  involved  in  thermal 
initiation  propagation,  the  reaction  model  can  be  greatly  simplified  with 
impunity.  In  Section  5  a  simple  model  is  presented  in  which  it  is  assumed  that  all 
of  the  energy  liberated  by  decomposition  in  the  initiated  region  of  the  material  is 
transported  immediately  to  the  front  and  used  up  there  to  sustain  the 
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propagation  of  the  front.  This  model  helps  fix  certain  useful  ideas  in  a  very 
simple  way. 

In  Section  6,  which  is  quite  lengthy  and  demanding,  the  fundamental 
hydrodynamic-reaction  equations  are  examined.  Comparison  with  experimental 
observations  allows  for  a  systematic  reduction  of  the  equations  until  finally  only 
a  thermal  diffusion-reaction  equation  remains.  This  is  the  basis  of  the  model 
equations  solved  in  Section  7  and  compared  with  experimental  results  in  Section 
8  for  the  purpose  of  estimating  model  parameter. 

The  establishment  of  the  thermal  diffusion-reaction  model,  which  was 
done  in  Section  6  is  believed  to  be  necessary  since  the  available  experimental  data 
is  too  sparse  to  validate  the  model  with  sufficient  conviction. 

The  analysis  shows  that,  with  sufficiently  detailed  experimental  data, 
reaction  parameters  above  the  initiation  threshold  can  be  obtained.  This  is  not 
generally  feasible  with  adiabatic  experiments  alone. 


a 


t  Experimental  Observation 


The  specific  physical  situations  to  be  modeled  are  the  reaction 
propagation  experiments  of  Stolovy,  Kidd,  and  Namenson.  Recent  experiments 
are  described  in  the  FY  90  report;  "Reaction  Propagation  Studies  in  Beam- 
Initiated  Confined  Explosives"[Ref  2].  Past  experiments  have  been  reanalyzed 
and  a  summary  is  given  in  an  appendix  to  this  report 

The  experiments  were  performed  on  short  cylinders  of  HE  (Figure  1) 
which  are  confined  in  an  aluminum  holder  (not  shown).  The  beam,  which  has 
been  collimated,  irradiates  only  the  central  region  of  the  HE  cylinder.  The  exit 
window,  by  which  the  beam 


Figure  1 

leaves  the  HE,  was  designed  to  be  the  weakest  part  of  the  holder  and  thus  the 
main  part  of  the  structure  to  be  disrupted  during  an  explosion. 

The  principal  objective  of  the  experiments  was  to  determine  if  the  HE  was 
thermally  initiated  beyond  the  beam  irradiated  region  and  to  measure,  by 
thermocouple  response,  the  time  of  arrival  of  the  initiation  front  at  several  radial 
positions.  Four  secondary  HE  were  studies:  TATB,  HBX-1,  PBXW-109,  PBX-9404. 

In  the  TATB  experiments,  explosions  were  weak,  with  no  initiation 
propagation  beyond  the  beam  region.  The  well  known  insensitivity  of  TATB  is 
evidently  retained  under  the  conditions  of  these  experiments.  In  the  HBX-1 
experiments,  the  samples  exploded  with  more  violence.  In  the  FY  90  set  of 
experiments,  no  initiation  beyond  the  beam  region  was  observed.  In  the  previous 
set  (see  Appendix  A.1)  in  which  the  beam  was  narrower  and  more  interne,  most 
(but  not  all)  of  the  HE  was  consumed;  a  propagation  velocity  of  the  order  of  1 
cm/sec  was  observed.  PBXW-109  behavior  was  generally  similar  to  HBX-1;  no 
quantitative  measure  of  the  propagation  velocity  was  obtained. 


Initiation  propagation  is  very  pronounced  in  PBX-9404  (94%  HMX).  In 
every  experiment  all  of  the  HE  was  consumed  and  a  violent  explosion  occurred. 
Initiation  propagation  velocities  (See  Appendix  A.1)  in  the  range  25-50  cm/sec 
were  observed.  Laser  experiments  (Appendix  A.1)  on  confined  PBX-9404,  in 
which  only  one  surface  of  the  sample  was  heated,  produced  similar  results. 
However,  surface  temperatures  greater  than  the  initiation  temperature  were 


necessary  to  precipitate  propagation. 

The  general  conclusions  that  have  been  drawn  from  the  analysis  are  as 
follows: 

1.  Initiation  propagation  is  an  intrinsic  property  of  confined  HE.  The 
values  of  the  frontal  velocities  vary  widely  from  one  HE  material  to  the  other. 
The  observed  representative  values  are:  TATB,  v  =  0;  HBX-1,  v  =  1  cm/sec;  PBX- 
9404,  v  m  25  cm/sec 

2.  In  those  HE  materials  in  which  propagation  is  present,  the  explosion  is 
delayed  until  after  significant  initiation  propagation  has  occurred.  This  indicates 
that  heat  and  probably  reactants  are  transported  out  of  the  initiated  region  into 
the  uninitiated  region.  Instead  of  creating  an  explosion  immediately,  the  energy 
goes  into  raising  the  temperature  of  the  adjacent  cold  region.  This  process  creates 
and  sustains  the  front  that  is  observed. 

3.  Confinement  is  necessary  for  thermal  initiation  in  general.  However, 
the  strength  of  the  confinement  does  not  appear  to  have  an  affect  on  initiation  or 
propagation.  This  indicates  that  the  function  of  confinement  is  the  retention  of 
gaseous  reactants.  The  strength  of  the  confinement  affects  the  violence  of  the 
eventual  explosion:  the  higher  the  pressure  at  which  explosive  disruption  occurs, 
the  greater  the  ensuing  violence. 

4.  The  evidence  is  consistent  with  the  concept  of  an  intrinsic  limiting 
propagation  radius  at  which  the  explosion  occurs.  In  TATB  this  radius  is  less 
than  0.6  cm.  In  HBX-1  and  PBXW-109  the  limiting  radius  is  probably  less  than  a 
centimeter.  In  PBX-9404  an  explosion  radius  of  about  1  cm  is  indicated. 

The  relatively  slow  propagation  speeds  and  weak  confinement  rule  out  a 
DDT  or  a  detonation  wave.  This  is  further  discussed  at  the  end  of  the  next 
section. 


2. 


2.  Conceptual  Physical  Description 


The  process  by  which  the  cylinder  of  HE  proceeds  from  beam  initiation  of 
its  central  region  to  eventual  explosion  can  be  conceptually  broken  up  into  the 
several  phases  depicted  in  Figures  2.  a-e.  This  construct  is  based  in  part  on 
experimental  observations  and  in  part  on  the  theoretical  and  modeling  analysis 
presented  in  Sections  6  and  7. 

Figure  2.a.  Initially  the  electron  beam  heats  the  cylindrical  volume 
bounded  by  radius  a.  The  temperature  Tj  is  less  than  the  initiation  temperature 
Tr.  Owing  to  oollimation,  the  beam  will  have  an  approximate  top  hat  profile. 
Normal  thermal  diffusion  produces  a  temperature  profile  with  an  rms  radius 
given  by 


(rJ)»2D0t  +  a*  (2.1) 

Normal  thermal  diffusivity  of  organic  HE  have  D0  <  10*3  an2/sec.  For  irradiation 
time  t  m  10  sec  and  oollimation  radius  0.2  cm,  the  temperature  profile  radius  will 
increase  by  <.04  cm  and  the  flat  portion  of  the  profile  will  decrease  by  the  same 
order  of  magnitude.  The  flat  portion  of  the  profile  will  reach  the  initiation 
temperature  in  an  adiabatic  manner.  In  modeling,  a  top  hat  profile  will  be  used 
with  an  effective  radius  determined  by  fitting  the  model  to  the  experimental 
data.  This  effective  radius  should  be  somewhat  less  than  the  collimation  radius. 

Figure  2,b,  Initiation  in  the  beam  region  has  been  achieved.  This  is  a 
transient  phase  in  which  the  conditions  for  propagation  are  developing. 
Necessary  conditions  are: 


a)  greatly  increased  thermal  diffusivity 

b)  creation  of  sufficient  heat  and  mobile  reactants  to  precipitate  frontal 
motion. 

The  temperature  rise  tends  to  be  adiabatic  before  propagation  begins.  The  time 
duration  of  this  phase  is  too  short  for  thermocouple  measurements.  In  PBX-9404, 
it  lasts  for  less  titan  a  millisecond .  If  propagation  does  not  develop,  a  thermal 
explosion  will  develop  on  this  time  scale. 

Figure  2  c.  Propagation  begins.  Heat  and  mobile  reactants  diffuse  from  the 
beam  region  to  the  front  Cold  material  is  heated  up  to  the  initiation  temperature 
Tr.  Below  Tr,  normal  diffusion  is  completely  negligible  on  a  millisecond  time 
scale.  Above  Tr  diffusion  is  very  pronounced.  The  temperature  in  the  central 
region  levels  off  as  a  result  of  the  greatly  enhanced  diffusion  and  the  consequent 
outflow  of  heat 


Figure  2  d.  A  quasi-steady  state  (QSS)  propagation  has  been  achieved,  the 
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temperature  profile  takes  on  a  Quasi-equilibrium  form.  Almost  all  of  the  energy 
generated  in  the  initiated  region  flows  to  the  frontal  boundary.  The  temperature 
profile  continually  adjusts  itself  to  the  motion  of  the  front  This  quasi-equilibrium 
status  is  a  consequence  of  the  high  thermal  diffusivity  that  develops  for  T>Tr. 

The  greatly  enhanced  thermal  conductivity  for  T  >  Tr  could  be  accounted 
for  by  the  following  hypothesis:  Once  the  initiation  temperature  Tr  is  reached  a 
profuse  generation  of  mobile  reactants  develops.  The  increased  mobility  accounts 
for  the  greater  thermal  conductivity.  The  mobility  of  the  reactants  (or  their 
diffusion)  is  curtailed  in  the  cold  region  beyond  the  front  Physical  changes  in  the 
initiated  region,  e.g.,  increased  porosity,  melting,  or  a  solid  state  phase  change, 
would  also  contribute  to  the  increased  diffusivity  in  this  region.  Although  there 
are  no  direct  experimental  observations  to  confirm  the  above  hypothesis,  such 
considerations  are  very  helpful  in  formulating  a  model  which  can  account  for  the 
observed  frontal  motion. 

Figure  2.e.  The  initiation  front  does  not  propagate  indefinitely.  As  the 
front  recedes  from  the  center  region  it  takes  an  increasingly  longer  time  for  heat 
generated  in  the  center  to  be  transported  to  the  front.  In  addition,  the  quasi- 
equilibrium  temperature  at  the  center  is  increasing  as  a  result  of  adjustment  to 
the  frontal  motion.  At  a  time  t*,  a  well-defined  frontal  position  R*  is  readied  for 
which  it  is  impossible  to  maintain  a  quasi-equilibrium  situation.  At  this  point  in 
time  the  temperature  in  the  central  region  rises  above  its  equilibrium  value  and, 
owing  to  the  autoacoelerating  nature  of  the  reaction  rate,  undergoes  an  explosive 
rise.  The  pressure  also  rises  explosively.  The  weak  confinement  experiments 
under  analysis  are  thus  abruptly  terminated  when  the  pressure  exceeds  the 
strength  of  the  confinement;  details  can’t  be  measured  by  thermocouple 
response. 

Additional  considerations  are  necessary  in  the  situation  (see  Figure  3) 
where  one  face  of  a  confined  HE  is  in  contact  with  an  inert  material  (usually  a 
metal)  which  is  at  a  temperature  greater  than  the  initiation  temperature  Tr  of  the 
HR  The  HE  is  initially  at  temperature  T0  <  Tr.  Experimentally  the  situation  has 
been  studies  for  PBX-9404  in  the  case  where  an  intense  laser  heated  the  inert 
material  (stainless  steel). 


Eventhough  the  interface  temperature  rises  above  Tr,  propagation  does 
not  begin  immediately.  Apparently  a  minimum  extent  of  the  HE  (indicated  in 
Figure  3  by  d)  needs  to  reach  the  temperature  Tr  before  initiation  propagation 
begins.  It  is  hypothesized  that  the  mobile  reactants  created  in  the  narrow 
interface  region  quickly  diffuse  into  the  cold  region  and  are  not  available  for 
further  reaction  at  the  interface.  When  the  point  d  attains  the  temperature  Tr 
initiation  begins.  Analysis  of  the  laser  experiments  for  PBX-9404  indicates  that 
this  HE  has  a  minimum  length  of  the  order  of  0.1  an;  the  particle  beam 
experiments  always  had  beam  radii  greater  than  this  and  thus  a  minimum  radius 
requirement  could  not  be  observed. 
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The  concept  present  here  and  the  analysis  to  follow  only  applies  up  to  the 
time  t*  when  the  front  reaches  the  explosion  radius  R*.  What  occurs  afterward 
does  so  on  a  time  scale  less  than  a  millisecond  and  depends  on  the  material.  The 
experiments  described  in  Appendix  A.1  indicate  that  for  TATB  no  further 
reaction  spreading  occurs;  for  HBX-1  additional  spreading  occurs  but  does  not 
consume  the  sample;  for  PBX-9404  additional  spreading  occurs  and  consumes 
the  whole  sample. 

For  the  purpose  of  distinguishing  thermal  initiation  propagation  from 
DDT  it  is  useful  to  describe  a  DDT  model  study  for  HMX  made  by  Mader  IRef  31. 
The  geometry  is  planar  and  the  extent  of  the  HE  is  greater  than  10  centimeters. 
The  confinement  is  very  strong  in  order  to  allow  a  deflagration  started  at  one  end 
of  the  HE  to  develop  into  a  detonation  wave  (c,  -  8  x  10= cm/sec)  at  the  other 
end.  In  Mader's  study  modeling  begins  when  several  centimeters  of  the 
deflagrating  end  of  die  HE  has  readied  kilobar  range  pressures;  a  detonation 
wave  develops  in  about  20  microseconds,  in  the  present  study  of  thermal 
initiation  propagation,  the  frontal  motion  is  measured  in  milliseconds  after 
initiation  is  first  achieved  in  the  beam  region.  Owing  to  the  relatively  weak 
confinement,  the  experiments  were  over  (i.e.,  disruption  occurred)  when  kilobar 
pressures  were  reached.  Since  the  propagation  time  scale  is  much  longer  than 
the  pressure  relaxation  time  (~p  secs),  isobaric  conditions  are  maintained;  large 
pressure  excursions  are  judge  only  to  occur  near  the  time  of  disruption. 


Figure  3 


Exothermic  Reaction  Modeling 


The  microscopic  chemical  reactions  which  govern  thermal  explosions  in 
condensed  energetic  materials  are  not  fully  known;  even  if  they  were,  they 
would  be  too  detailed  for  performance  applications.  In  practice,  calculations  are 
often  based  on  model  reactions  which  utilize,  as  variables,  the  fraction  NAof 
molecules  of  the  material  which  remains  unreacted,  and  Nb  the  fraction  that  has 
reacted.  The  best  models  are  based  on  actual  experiments  and  have  their  limits  of 
applicability  defined. 

The  simplest  model  corresponds  to  unimolecular  decomposition.  The 
reaction  scheme  is 


A^B  +  A';Q  (3.1) 

The  molecule  A  decomposes  thermally,  with  rate  constant  k,  into  two  (or  more) 
different  molecules  and  an  amount  of  heat  Q  (per  unit  mass)  is  released.  For  a 
homogeneously  reacting  system  the  corresponding  rate  equation  is 


^  =  -kA(T)NA  (3.2) 

The  temperature  dependence  of  the  *ate  constant  usually  (though  not  always) 
takes  the  Arrhenius  iorm: 


KA(T)=*Zexp(-TA/T)  (33) 

The  pre-exponential  factor  Z  has  units  of  inverse  seconds  and  T  is  the  absolute 
temperature.  The  activation  temperature  is  defined  as  follows: 


TA  =  EA/Rg  (3.4) 

where  EA  is  the  activation  er  orgy  (cal/mol)  and  Rg  is  the  gas*  constant  Rg  ■  1.987 
cal/mol-deg.  The  fraction  of  unreacted  molecules  is  defined  as 

=  (3,5) 

where  n©  is  the  original  number  of  molecules  per  unit  volume  and  nA  is  the 
number  of  molecules  which  are  unreacted.  The  rate  at  which  exothermic  energy 
is  released  is 

45  — (3.6) 
at  at 

in  which  q  is  the  heat  per  unit  mass  acquired  by  the  material. 
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A  more  detailed  model  corresponds  to  an  autocatalytic  scheme: 


A  b  B  +  A';  Q 

(3.7a) 

B  +  A  bB+B+A;  Q 

(3.7b) 

The  molecule  A  decomposes  into  a  molecule  B  and  other  products;  then  the 
molecule  B  reacts  with  other  molecules  A  to  produce  a  decomposition  similar  to 
the  previous  one,  without  itself  changing.  The  first  reaction  occurs  thermally,  i.e., 
as  a  result  of  elevated  temperature;  the  second  is  similar  but  occurs  because  of 
the  presence  of  B  which  acts  cataly  tically,  i.e.,  it  is  not  changed  in  the  reaction. 
The  corresponding  rate  equations  are: 


^  =  -kA(T)NA-kB(T)NANB  (3.8  a) 

^-MT)NA+kB(T)NANB  (3.8  b) 


in  which 


N„»S!  (3.9) 

and  ne  is  the  number  of  reacted  molecules.  It  is  convenient  to  introduce  the 
notation 


KlNjbNBJTi-k^TJNA  +  kBClINANB  .  (3.10) 

The  rate  of  heat  generated  per  unit  mass  decomposed  is 

^-Q^-QKlN^NBiT)  .  (3.11) 

The  reaction  Eqs.  (3.8  a,  b)  imply  the  following  conservation  relationship: 

Na  +  Nb®  const  (3.12) 

If  at  the  start  of  an  experiment  (t  =  0)  only  type  A  is  present  then 

NaI^O)®!  ,  NB(t*0)*0  ,  (3.13) 


11 


thus 


Na  =  1-Nb 


(3.14) 


and 


K  =  kA(l-NB)+kB(l-NB)NB  #  (3.i5) 


Rogers,  Janney,  and  co-workers  [Ref.  4]  have  studied  a  large  number  of 
explosives  by  isothermal  differential  scanning  calorimetry  (DSC)  and  have 
analyzed  their  results  with  a  generalized  autocatalytic  reaction  scheme: 


^ = k  a(X-Nb)" + kB  (1-  NbJPNb’ 

They  find  that  the  exponents  have  the  value  unity,  at  least  up  to  50% 
decomposition,  and  in  some  cases  beyond.  Temperature  variations  yield  reaction 
rates  kACI)  and  kBCD  that  have  the  Arrhenius  form  over  a  wide  temperature 
range;  however  the  constants  differ  for  each  well  defined  region. 

The  temperature  ranges  studied  by  the  DSC  method  except  for  TATB,  do 
not  go  up  to  the  thermal  initiation  temperature  region  investigated  by  Stolovy,  et 
al.  This  is  simply  because  the  exothermic  reaction  rates  become  too  large  to  hold 
the  DSC  system  at  constant  temperature. 
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4.  Adiabatic  Thermal  Initiation  by  Partide  Beams 


The  thermal  initiation  of  energetic  materials  by  partide  beams  (electrons 
and  protons)  has  been  studied  [Ref  1.]  primarily  in  an  adiabatic,  homogeneous 
context  with  confinement.  The  material  is  uniformly  irradiated  by  a  partide 
beam  which  deposits  energy  uniformly  at  the  rate  qb  per  unit  mass.  The  pertinent 
equations  are 


^4=-k(na^jb;T) 

(4.1a) 

^  =  K(Na^b;T) 

(4.1  b) 

|i-Cv|t  =  4b  +  QK(NAJJB;T) 

(4.1c) 

The  equations  imply 

d(NA  +  NB)  a 
dt 


(42a) 


(42b) 


integrating  from  room  temperature  (t  ■  0,  T  «  To,  NA(0))  yields 


Na-i-Nb 

(43a) 

fCvdT-QNB-qbt 

(43b) 

The  last  equation  simply  states  that  difference  between  the  internal  energy  and 
the  decomposition  energy  is  supplied  by  the  beam  energy  deposition. 


Experimentally,  the  temperature  as  a  function  of  irradiation  time  is  the 
primary  data.  Generic  behavior  is  depicted  in  Figure  4: 
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The  initiation  temperature  is  dearly  indicated  by  the  radical  slope  change  at 
T=Tr.  Below  Tr,  evidence  for  exothermic  reactions  is  not  very  pronounced,  while 
above  Tr  it  is  overwhelming.  The  amount  of  deposited  beam  energy  required  to 
bring  the  HE  to  initiation  is  called  the  initiation  energy  (from  25°  C)  and  is 
denoted  by  Er.  Clearly 

ER-qbtR  (from  25°  C)  (4,4) 

Extensive  measurements  of  Tr  and  Er  have  been  made  by  Stolovy,  et  al  [Ref  1], 
For  PBX-9404,  whose  propagation  properties  will  be  analyzed  quantitatively  in  a 
later  section,  the  values  are 

Tr*  28015°  C  (553  ±5°K) 

(45) 

Er  *  635  ±  05  cal/gm  (from  25°  C) 

The  beam  heating  rate  in  these  experiments  varied  between  7.62  and  11.28 
cal/gm-sec,  with  an  average  at  about  10  cal/gm-sec. 

The  detailed  temperature  behavior  in  the  region  around  Tr,  when  looked 
at  under  high  resolution,  shows  considerable  variation  from  one  HE  material  to 
another.  For  PBX-9404  the  temperature  achieves  a  plateau  just  before  initiation; 
this  indicates  a  phase  change  which  is  consistent  with  the  known  melting  at 
2S5°C  TATB  exhibits  a  smooth,  but  rapid  change  at  Tr  »  409°C.  The  melt  occurs 
at  Tm  -  450°C.  The  behavior  of  HBX-1,  which  consist  of  40%  RDX,  35%  TNT  and 
~16%  AX,  is  quite  complex.  Examples  are  given  in  Ref.  1. 

Attempts  to  analytically  understand  thermal  initiation  by  particle  beams 
through  use  of  reaction  models  with  parameters  taken  from  relatively  tow 
temperature  thermal  decomposition  experiments  have  not  been  successful  in 
reproducing  the  observed  details  in  the  vidnity  above  and  below  Tr.  This  would 
be  explained  if  the  reaction  rate  underwent  a  radical  change  at  Tr.  In  this  case  the 
function  of  the  partide  beam  is  to  (uniformly)  bring  the  HE  up  to  die 
temperature  region  where  the  radical  change  occurs.  The  initiation  temperature 
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then  is  an  intrinsic  property  of  the  HE.  In  what  follows  it  will  be  seen,  at  least  for 
PBX-9404,  that  the  initiation  front  propagates  with  a  temperature  which  is  the 
same  as  the  measured  adiabatic  initiation  temperature,  however  the  associated 
exothermic  rate  is  several  orders  of  magnitude  higher  than  the  beam  rate  (  qb~10 
cal/gm-sec)  used  in  the  adiabatic  experiments.  This  lends  support  to  the  concept 
of  a  radical  change  occurring  at  a  temperature  Tr  which  is  intrinsic  to  the 
materials. 

For  analysis  of  propagation  the  reaction  rate  needs  to  be  known  in  a 
temperature  region  <  40°C  above  the  initiation  temperature  Tr.  The  following 
form  will  be  used 


KQsqRexpl^T-TR), 


(4.6) 


in  which  qR  is  the  exothermic  rate  at  initiation  and  Ta  is  an  effective  activation 
temperature.  In  Appendix  A.2  it  is  shown  that  the  autocatalytic  reaction  rate  Eq. 
3.10,  with  kA(T)  and  ka(T)  in  the  Arrhenius  form,  assume  the  form  in  Eq  (4.6)  for 
a  limited  temperature  range.  For  given  reaction  constants  in  Eq.  3.10,  Za,  Zb,  Ta, 
Tb,  the  quantities  qR  and  To  can  be  determined;  lacking  the  former,  the  latter 
quantities  must  be  treated  as  model  parameter  to  be  determined  by  comparison 
with  experiment.  The  value  of  the  initiation  energy  and  temperature  from 
adiabatic  experiments  are  important  parameters  entering  into  the  propagation 
equation. 


5.  Simple  Model  for  Initiation  Propagation 


In  order  to  fix  notions  for  the  dynamics  of  the  moving  front  it  is  useful  to 
construct  a  simple  model:  it  is  assumed  that  diffusion  is  small  for  T  <  Tr  and  that 
heat  transport,  by  any  means,  is  very  large  for  T  >  Tr.  Once  initiation 
propagation  begins,  the  temperature  in  die  initiated  region  remains  constant 
since  all  of  the  newly  generated  energy  moves  immediately  to  the  cold  region. 
Figure  5  illustrates  the  model  situation. 


Tit. 


Tc 


3.  R 

Figures 


As  a  matter  of  convenience,  the  beam  is  turned  off  once  initiation  is  achieved. 
The  energy  increment  of  the  initiated  region  is 


U-U0»ErMr  (5.1) 

Where  Mr  is  the  mass  of  the  initiated  region  and  Er  is  the  energy  per  unit  mass 
required  to  take  the  material  from  T0  to  Tr.  The  experimentally  measured 
initiation  energy  will  be  used  in  calculations.  The  rate  at  which  energy  changes  in 
the  initiated  region  is  given  by 


^-SrMr  (5.2) 

in  which  qR  is  the  rate  of  change  of  energy  per  unit  mass.  This  quantity  is  not 
assumed  to  be  the  same  as  the  beam  energy  rate  which  produced  initiation.  Both 
qR  and  Tr  are  quantities  established  when  initiation  spreading  occurs. 

Energy  conservation  requires 

^^rMrsEr^^  (53) 
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For  cylindrical  geometry,  with  R  =  radius  and  1  =  length. 


M  ss  eR2  1  |A 

where  p.  is  the  mass  density.  Equation  (53)  reduces  to 


(5.4) 


c1r_  2  dR 

Er  '  R  dt 


(53) 


The  time  constant  tj  =  ER/qR  is  the  tine  required  to  adiabatically  change  the 
energy  of  a  unit  mass  of  the  material  by  an  amount  Er  if  the  heating  rate  isqR. 
The  solution  to  Equation  (53)  is 


1 

ti 


(5.6) 


in  which  it  is  assumed  that  propagation  begins  at  point  a  at  time  t  =  0.  This 
relation  will  now  be  applied  to  the  PBX-9404  data  found  in  the  appendix  and 
repeated  here: 

acou  *  0.238  cm  (collimation  radius) 

Ri  ■  0.476  cm  tj  *  25.0  msec  (5.7) 

R2  *  0.953  cm  tj  »  433  msec 

Assuming  that  propagation  begins  at  the  ooilimation  radius,  the  data  at  position 
2gives 


,i*2E^srls'8m“c  <s® 

Alternatively,  the  use  of  the  information  at  the  two  positions  gives 

t, «  U-JA  ^  *  133  msec 

In  a  *  In  R2-^* -1.67  (5.9) 

a  =  0.188  cm  (model  value) 

The  values  obtained  have  about  a  20%  uncertainty. 
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If  now  the  experimental  value  Er  =  64  cal/gm  is  used,  then 


qR--£“= 4.74  x  103  cal/gm-sec  (5.10) 

The  large  magnitude  of  this  value  relative  to  the  beam  heating  rate  (~  10  cal/gm* 
sec)  is  not  changed  by  subsequent  refinements  in  the  analysis.  It  is  simply  the 
minimum  rate  at  which  exothermic  energy  must  be  generated  in  order  to 
produce  a  front  with  the  observed  motion. 

The  data  in  Equation  (5.7)  indicate  that  the  explosion  time  after  initiation 
has  occurred  is  something  in  excess  of  44  msec  It  is  instructive  to  compare  this 
to  the  explosion  time  in  the  adiabatic  case,  i.e.,  no  propagation.  A  simplified 
exothermic  heating  rate  model  equation  is  (see  Section  4 ): 

Cv  **  Tr)  (5.11) 

The  time  to  explosion,  i.e.,  the  time  it  takes  for  the  temperature  to  go  from  Tr  to 
oo,  is  given  by 

*.4*f  expf-^(T-TR)|dT  =  ^Si  (5.12) 

q Rj*  l  V  qR  »a 

For  PBX-9404,  a  standard  value  of  Ttt  and  Cv  are 

Ta*  26,000*  K 


Cv»  0.4  cal/gm  (high  temp.)  (5.13) 

These  give,  with  the  experimental  Tr  «  553°  K, 

t,  *  1  msec  (5.14) 

Thus  the  actual  explosion  time,  after  initiation,  is  increased  by  more  than  an 
order  of  magnitude  as  a  result  of  propagation. 

It  is  also  instructive  to  calculate  the  fraction  of  the  original  material  which 
must  decompose  in  order  to  generate  the  energy  that  keeps  the  front  moving,  this 
is  given  by 


amb!e 

M  Q 


(5.15) 


in  which  Q  is  the  exothermic  energy  per  unit  mass  decomposed  and  Er,  the 
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initiation  energy,  is  the  energy  required  to  bring  the  cold  material  at  the  front  to 
die  initiation  temperature  Tr.  The  extreme  values  of  Q,  for  PBX-9404,  are  given  as 

_  /  500  cal/gm  (Thermal  Explosion)  \ 

'*~\1500cal/gm  (Detonation)/  (5.16) 

Since  Er  =  64  cal/gm,  the  above  yields 


(5.17) 

Thus,  the  order  of  10%  of  the  available  exothermic  energy  is  expended  on 
propagation;  the  remainder  will  go  into  the  final  explosion  when  it  occurs. 

The  main  failing  in  the  simple  model  considered  here  is  that  propagation 
is  unlimited.  Inversion  of  Equation  (5.6)  yields 

R(t) «  a  exp(t/2  tj)  (5.18) 

The  more  elaborate  model  to  follow,  in  which  finite  diffusivity  and  mass  flow  is 
accounted  for,  contains  naturally  a  limiting  radius  R*  at  which  an  explosion  will 
occur.  For  finite  diffusivity  there  will  always  be  a  point  at  which  all  the  energy 
generated  in  the  vicinity  of  R  »  0,  cannot  diffuse  to  the  front.  The  temperature  in 
the  central  region  will  then  rise  and,  owing  to  the  autoaccelerating  nature  of  the 
reactions  rate,  will  quickly  destablize  the  quasi-equilibrium  state  essential  for 
propagation.  For  PBX-9404,  the  value  of  R*  will  be  a  centimeter.  However  the 
value  of  the  diffusivity  will  be  -103  cm2/sec.  Owing  to  this  very  large  diffusivity, 
the  values  of  tj,  a,  and  %  obtained  in  the  simple  model  are  very  nearly  correct; 
only  the  extent  of  propagation  is  modified.  The  extent  of  propagation,  R*,  is  very 
important  since  it  determines  the  minimum  amount  of  HE  which  will  undergo 
thermal  explosion. 


•MS^SO. 
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In  this  section  general  transport  equations  and  concepts  will  be  developed. 
The  general  analysis  techniques  are  based  on  Landau  and  Lipchitz's  "Fluid 
Mechanics"[Ref.  5].  Thermal  initiation  propagation  modeling,  based  on 
experimental  observations,  will  be  then  performed. 

The  basic  conservation/generation  equation  has  the  general  form: 

|£  +  V-)0  =  S$  (6.1) 

in  which  0  is  a  density  (e.g.,  mass,  momentum  component,  energy,  etc,  per  unit 
volume);  is  the  corresponding  flux  vector.  If  Spanishes,  and  if  vanishes  at 
the  boundary  of  the  space  under  consideration,  then 


J  $  d  u  =  const  (62) 

i.e.  the  total  amount  of  the  quantity  in  the  volume  u  is  conserved.  When  a 
quantity  is  not  conserved  its  rate  of  generation  per  unit  volume  is  specified  by 
the  source  term  S*. 

The  HE  material  is  assumed  porous.  When  die  initiation  temperature,  Tr, 
is  readied  the  material  decomposes  exothermically  and  produces  solid  and 
gaseous  products.  Physical  changes  in  the  solid  substrate  may  precipitate 
initiation.  Hot  gases,  filling  the  pores,  produce  increasing  pressure  which,  if  not 
relieved,  will  produce  an  explosion  in  less  than  a  millisecond.  Relief  can  be 
obtained  by  heat  conduction  and  mass  flow  out  of  the  initiated  region  into  the 
cold  uninitiated  region.  Since  pressure  relaxation  takes  place  on  a  time  scale 
much  shorter  than  the  observed  initiation  propagation  time,  pressure  wave 
propagation  will  not  be  significant  until  the  explosion  occurs.  Mass  transport  of 
the  gaseous  products  through  the  porous  materia!  will  be  allowed  for.  General 
thermal  conduction  is  taken  Into  account  and  will  be  shown  to  dominate. 


The  thermal  initiation  propagation  model  to  be  developed  is  meant  to 
apply  from  the  time  when  the  initiation  temperature  is  reached  until  the  time 
when  the  quasi-steady  state  propagation  can  no  longer  be  maintained.  At  this 
latter  time  rapid  temperature  and  pressure  changes  begin  and  an  explosion 
develops  on  a  millisecond  time  scale. 

Specific  reaction  modeling  will  be  a  generalization  of  the  autocataiytic 
reaction  process  to  include  energy  and  mass  transport  For  convenience  the 
process  is  repeated  here: 
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A  ->  B  +  A* 

B  +  A->2B  +  A' 


(63) 


The  product  B  is  assumed  to  be  a  gas  and  A'  is  assumed  to  be  a  solid  (or  at  least 
an  immobile  liquid). 

A  microscopic  description  of  the  phenomena  to  be  modeled  is  not 
available.  It  is  thus  appropriate  to  take  as  the  elementary  material  "point"  a  small 
macroscopic  volume  which  contains  many  pores.  The  motion  of  the  gas  within 
the  elementary  "point"  is  locally  convective  and  highly  chaotic.  No  attempt  will 
be  made  to  model  behavior  within  the  "point".  The  average  properties  over  the 
small  macroscopic  volume  will  be  used.  These  will  either  be  inferred  from 
experiment  or  postulated  by  reasonable  argument 

For  each  component  i  *  A,  A‘,  B  the  number  per  unit  volume  is  designated 
by  nu  the  initial  value  of  nA  is  no.  It  continues  to  be  convenient  to  express 
numbers  of  molecules  as  fractions  of  the  initial  number  of  type  A: 

=  (6,4) 

The  two  solid  components,  A  and  A*,  and  their  molecular  masses  are  related  as 
follows: 


nA*n0-nA 


(63) 


Through  these  relations  any  explicit  dependence  on  A' can  be  removed.  It  should 
be  noted  that 


nA^na  (63) 

since  8  is  free  to  move  from  its  plats  of  origin,  whereas  A'  is  not 


Although  number  density  is  most  convenient  in  discussing  reactions,  mass 
density  must  be  introduced  for  hydrodynamics.  The  following  relations,  for  mass 
density  m,  hold: 


Pi «  mi  n|  s  PoN* 

where 


(67) 


Po«mAiio 


(63) 


21 


The  total  solid  density  is 


(6.9) 


and 


d|is  =  ^HodN*  (6.10) 

The  fractional  volume  occupied  by  the  pores  is  defined  as 

(6.11) 

where  v  is  the  total  volume.  The  quantity  is  assumed  not  to  change  significantly 
in  the  initiated  region  during  the  quasi-steady  state  propagation  period.  The 

average  macroscopic  gas  density  4b  is  related  to  the  true  average  gas  density,  pe , 
in  the  pores  by 


4b»cpb  (6.12) 

The  fundamental  kinematic  variable  is  the  average  mass  fiow  velocity  V 
which  is  defined  as  the  average  gas  momentum  per  unit  gas  mass.  The 
momentum  per  unit  macroscopic  volume  is  then  given  by  4b^*  The  mass 
conservation/generation  equation  for  Bis  then 

+v.ne\?.p<£  (6.13) 


in  which  pJR  is  the  rate  of  generation  (per  unit  macroscopic  volume)  of  gaseous 
mass  which  results  from  reactions  in  or  with  the  solid  substrate.  It  follows  that 


3Nb 

at 


+ 


V-Nd? 


In  Section  3  the  reaction  rate  for  Nb  was  designated  by  K,  thus 


For  the  solid  components,  neither  of  which  move. 


3N v  8nA 


«-K 


(6.14) 


(6.15) 


(6.16) 
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and 


=  HoR  (6.17) 

Turning  now  to  momentum  transfer#  which  by  previous  assumption  (i.e., 
no  pressure  waves  in  the  solid)  will  only  involve  the  gas#  the 
conservation/generation  equation  is  (neglecting  internal  gas  viscosity) 

+  (pBV\UPl}*pBf  (6.18) 

in  which  P  is  a  pressure  and  1  is  the  unit  tensor.  This  equation  is  first  obtained 
in  terms  of  the  true  density  and  pressure  in  the  pores#  pe  and  P,  and  then 
multiplied  by  o.  Thus 


P  =  crP  (6.19) 

The  only  source  (actually  a  sink)  of  gas  momentum  is  the  reduction  of  the  flow 
momentum  through  interaction  with  the  porous  substrate.  The  force  per  unit 

mass#  ?  #  results  from  viscous  effects  and  flow  randomization  that  occur  within 
the  material  "point";  it  is  reasonably  (and  usually)  modeled  by 

(620) 


in  which  Q#  the  flow  resistance,  will  be  assumed  constant  during  the  propagation 
period. 

Eq.  (6.13)  implies  the  following  identity  for  an  arbitrary  quantity  G: 

— g--—  +  V  -(ligG  +  (621) 

in  which 

jjL»|.  +  9.V  (622) 

at  at 

Setting  G  *  9  in  Eq.  (621)  and  combining  the  result  with  Eq.  (6.18)  results  in 

Pe^s-VP+pel-poKT^  (623) 

Although  the  generation  of  gas  (which  has  2ero  average  velocity)  is  not  a  source 
of  momentum#  the  newly  generated  gas  acquires  momentum  from  the  flowing 
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gas;  this  causes  a  reduction  in  the  average  velocity  which  is  accounted  for  by  the 
last  term  in  Eq.  (6.23). 

In  making  energy  consideration  *,  the  system  as  a  whole  is  taken  as  closed, 
i.e.,  no  energy  is  introduced  from  the  outside;  the  beam  has  either  been  turned 
off,  or  ( as  was  seen  in  Section  5)  the  beam  energy  deposition  rate  is  negligible 
compared  to  the  reaction  energy  rate.  The  conservation  equation  for  the  system 
as  a  whole,  i.e.,  solid  plus  gas,  is 

^  4  V '  Mh  *  +  ^  H  =  0  (6.24) 

in  which  es  and  eg  are  the  solid  and  gaseous  energies  per  unit  mass,  and 

hB“eB  +  ]I§  =68  + (6.25) 

is  the  enthalpy  per  unit  mass  of  the  gas.  The  heat  flux  vector «}  accounts  for  any 
energy  transport  that  is  nu  associated  with  mass  transport. 

Several  transformation  will  now  be  Hade  which  will  make  internal 
heating  effects  (and  entropy  transport)  explicit.  First  apply  Eq.  (6.21)  with  G  = 

hB  +  ^  V2.  When  combined  with  Eq.  (6.24)  this  results  in 

^  +  ^-¥  +  W^2'?1  +  ,1<'E(hB  +  292)tV ^  =  0  <M6> 
Next ,  take  the  scalar  product  of  V  with  eq.  (6.23).  The  result  is 

«-y.VP+jiBV-f-*ioK$2  <&27) 

Inserting  this  into  Eq.  (6.26)  yields 

i4hB+^l)+v.q=o  (628) 

Utilization  will  now  be  made  of  the  following  thermodynamic  relations  which 
hold  for  either  the  solid  or  the  gaseous  component: 

dEi®Tdsi  +  (6.29a) 
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(6.29b) 


dh,  =  Tdsi  +  ^ 


d(m£i)  =  |ii  T  d  Sj  +  hj  djii 


(629c) 


in  which  S\  is  the  entropy  per  unit  mass  of  the  Ith  component.  Introducing  these 
into  Eq.  (628)  results  in 


^^  +  WT(^  +  ?VsB]^0R(hs-hB-l92) 

-iib^  (?-^R^)-v.  q=o 


(6.30) 


This  is  the  general  equation  of  heat  , transfer  for  the  solid-gas  system  undergoing 
exothermic  reactions.  The  first  term  on  the  left  is  the  rate  of  heating  (per  unit 
volume)  of  the  solid;  the  second  term  is  the  rate  of  heating  of  the  moving  gas.  The 
right  hand  side  exhibits  the  internal  sources  of  heat  The  first  term  results  from 
decomposition  with  the  non-heating  kinetic  energy  subtracted  off;  the  second 
term  is  the  heating  rate  produced  by  the  dissipation  of  kinetic  energy  of  the 
flowing  gas  as  a  result  of  the  interaction  with  the  porous  solid  substrate  and  the 
wwly  created  (zero velocity)  gas;  the  last  term  arises  from  thermal  conduction, 
tue  decomposition  term  can  be  written  in  terms  of  the  heat  of  reaction  as  follows: 

^(h.-hs-l^J.HoKfQ-^V2)  (631) 


where 

Q«af(h,-h„j  (632) 

Seat  flux  vector,  on  general  grounds,  is  taken  to  be 

q«-icVT  (633) 

The  thermal  conductivity  tc  is  assumed  not  to  vary  significantly  for  the  quasi¬ 
steady  state  variation  to  be  considered  (A  T  $  40°  C )  .  The  heating  rates,  as 
usual,  are  written 


8$i  3siat  8T 

T  at  ”t3t  aT*'1  at 


(634) 


in  which  cj  is  an  appropriate  specific  heat  A  detailed  analysis  is  not  in  order;  it 
will  suffice  to  write 
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(6 35) 


dss  dse  8T 

^sr^ar^aF 


where 


CyS 


M>  sCs  +  ^BCB 

m  +  HB 


(636) 


is  an  effective  specific  heat  and 


+  (6 37) 

is  the  total  mass  density.  Whenever  a  value  for  Cv  is  needed,  an  estimated  high 
temperature  value  for  the  unreacted  HE  (the  dominant  component)  will  be  used. 
Eq.  (630)  can  now  be  written 

*v¥+^t9Vsb^k(q-^v>) 

-HB^-p-]^KV)-v.q=o 

Having  now  established  sufficiently  general  equations  governing  mass, 
energy,  and  momentum  transfer  for  the  system,  the  analysis  will  now  turn  to 
comparing  the  roles  of  mass  flow  and  heat  conduction  in  transporting  energy 
toward  the  cold  region  during  the  quasi-steady  state  thermal  initiation 
propagation.  Of  necessity,  the  analysis  is  done  specifically  for  PBX-9404,  however 
general  conditions  are  stated  so  that  applications  will  be  facilitated  whenever 
sufficient  information  is  available  for  other  HE  materials. 

the  main  factors  to  be  used  are: 

1.  The  flow  velocity  must  not  be  greater  than  the  initiation  front 
velocity,  then 


V<  SO  cm/sec  (639a) 

1/2  V*  <  135  x  10*  ergs/gm  (639b) 

This  factor  has  already  been  used  to  rule  out  sound  (i.e.  pressure)  propagation. 
The  solid  components  are  considered  mechanically  rigid. 

2.  the  heat  of  reaction  is  quite  large: 

Q -500 -1500  caJ/gm 
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-  2  x  1010  -  6  x  1010  ergs/ gm 


(6.40) 


Since,  in  the  sample  model  of  Section  5,  the  rate  of  heat  production  was  deduced 
to  be  quite  high,  the  larger  possible  values  of  Q  are  favored.  A  value  of  Q  *  1000 
cal/gm  =  4  x  1010  ergs/gm  will  be  used  in  estimates 

3.  The  rate  of  increase  of  the  heat  in  the  initiated  region  during  the 
propagation  period  must  be  considerably  less  than  the  exothermic  reaction 
energy  rate: 


\l  «^oQK  (6.41a) 

If  this  did  not  hold,  an  autoaccelerating  temperature  increase  would  bring  about 
an  explosion. 

A  fundamental  conclusion  is  arrived  at  by  looking  at  Eq.  (638)  in  the 
vicinity  of  the  center  of  cylinder  (r  =  0).  Since  the  mass  flux  must  vanish  here,  the 
velocity  also  vanishes.  Thus  Eq.  (6.38)  reduces  to 

HCv~«-V§»HoKQ  (6.41b) 

Eq.  (6.41a)  then  implies  that  thermal  conduction  dominates  the  region  about  the 
origin.  Using  Eq.  (633),  Eq.  (6.41b)  can  be  written 

~-DtV2T»^P§;K  (6.41c) 


where 


Dt  *  (6.41  d) 

is  the  thermal  diffusivity  (units  -  cmVsec).  Once  initiation  is  achieved  in  the 
beam  region,  an  explosion  can  be  prevented  by  thermal  conduction  if  the  heat 
diffuses  out  of  the  beam  region  (r  £  a)  on  the  time  scale  t»  =  1  msec  for  an 
adiabatic  explosion.  This  gives 


>a 


or 


(6.41e) 


Dir>  .25  x  10*  cm2/sec 
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Thus  a  very  large  thermal  diffusivity  is  anticipated. 


In  due  course  all  velocity  dependent  terms  in  Eq.  (638)  will  be  analyzed. 
To  begin  with,  Eqs.  (6.39b)  and  (6.40)  imply. 

V*S|v2<<Q  (6.42) 

Thus  Eq.  (638)  becomes 
3T 

UC^^  +  jlBTV.VsB  +  V.qs^KQ-liB^f  (6.43) 

Further  reduction  will  come  after  consideration  of  the  momentum  transport. 

For  one-dimensional  flow,  the  momentum  flux  contained  in  Eq.  (6.18)  is 

jnuwn  “  M-B  (v2  +  jig)  '  (6.44) 

The  gas  will  be  assumed  tc  be  ideal  Thus 

P«anekT  (6.45a) 

If  the  product  B  is  a  single  molecule  then  a  *  1.  If  the  product  consist  of  a 
molecules ,  their  partial  pressure  add  up  to  give  Eq.  (6.45a).  It  is  assumed  that  the 
relative  concentrations  do  not  change  during  propagation.  Introducing  mass 
densities  yields. 

**  ^  ^  (6.45b) 

In  which  Rg  *  83  x  107  ergs./mokieg  (gas  constant).  Since  T  >  553°  (=Tr)  and  the 
average  molecular  weight  Mb/ a  <  300  gms/mol  (HMX),  the  following  holds: 

£  >1.5x10*  cmVsec2  >  >  V2  (6.46) 

Using  this  result  and  Eq.  (6.20),  Eq.  (6.18)  becomes 

V  - 

— gp  +  (6.47) 
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This  equation  is  formally  solved  as  follows: 


|Ab^  =  M>b^  J  efl(t'-t,)  VPdt'  (6.48a) 


Repeated  integration  by  parts  yields 


=  /  ,_e-  Q  -  v(  £  -  |£  h-  O  (^ )]  (6,48b) 

For  a  sufficiently  large  flow  resistance  ft.  Eq.  (6.48)  reduces  to 

=  (6.49) 


Thus,  in  a  time  interval 


At j~Q“*  (6.50) 

The  mass  flux  relaxes  to  a  quasi-steady  state  value  given  by  Eq.  (6.49).  Once  the 
rate  of  pressure  increase  becomes  explosively  large  this  situation  is  terminated. 

The  flow  resistance  ft  is  the  primary  parameter  determining  the  mass 
flow.  Unfortunately,  the  experiments  do  not  give  any  measure  of  mass  flow.  The 
frontal  velocity,  which  marks  the  progress  of  initiation,  gives  only  a  possible 

upper  bound  Despite  this  uncertainty,  a  useful  lower  bound  to  ft  can  be 
conjectured  as  follows:  It  is  reasonable  to  assume  that  the  drag  force  due  to  the 
porous  material  is  much  greater  than  the  drag  force  due  to  the  creation  of  zero 
(average)  velocity  gas. 

Thus,  from  Eq.  (6.23), 


(651a) 


or 


Q>>& 


From  the  simple  model  presented  in  Section  5, 


K 


Sr 

Q 


(651b) 


(652a) 


29 


Thus 


AM  =  Nb,NbQ=Er 


(6.52b) 


q»SR=1 
Er  tj 


(653a) 


or 


Q“ 1  s  A  tf « ti  *  12  msec  (653b) 

Thus  a  flow  relaxation  time  of  the  order  of  a  millisecond  is  not  unreasonable  and 
assures  that  the  quasi-steady  form  given  by  Eq.  (6.49)  holds  during  the 
propagation  period  which  lasts  about  40  msecs. 


Eq.  (6.49)  now  replace  the  momentum  equation.  The  mass  and  energy 
equations  now  take  the  form 


H  <5  «■  HB  $ •  (T  Vsb  ♦  V  •  q  -  PoQR 


(654) 

(655) 


The  effort  will  now  be  directed  toward  deriving  conditions  under  which 
the  velocity  dependent  term  in  Eq.  (6.55)  can  be  neglected  relative  to  the 
conduction  term.  Applying  Eq.  (&29b)  to  the  present  situation. 

TVSB  +  ^»VhB  (656) 


For  an  ideal  gas,  in  a  limited  temperature  range, 


thus, 


dhsCpdT 

+  CpeV  T  +  V  q-PoKQ 


(657) 

(658) 


Reintroducing  cj « -  k  V  T,  and  redefining  a  thermal  diffusivity 


Eq.  (658)  becomes 
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“'*1 +V.g  =  HoK<2 


(638) 


J 


Appropriately  neglecting  the  time  derivative,  a  formal  solution  of  eq.  (6.58),  for 
cj,is(in  cylindrical  coordinates): 


r  q  =  ev'£  e"vr'  \i0  K  Q  r  d  r*  (6.59) 

where 


v  = 


iB  cpB  V 

£  cv  Dt 


(6.60) 


If  the  exponential  in  Eq.  (6.59)  are  essentially  unity,  the  velocity  dependent  term 
in  Eq.  (6.58)  is  negligible.  This  requires 


vr«l 


(6.61a) 


or 


(6.61b) 


The  maximum  r  is  about  1  cm.,  the  maximum  v~50  cm/sec.  The  following  is 
considered  reasonable. 


•Ol<0T^<1  (6.62) 

Conservatively  then, 

Dr>>50cm2/sec  (6.63) 

is  a  sufficient  condition  for  neglecting  the  velocity  dependent  term  in  Eq.  (6.58). 
Comparing  this  with  Eq.  (6.41e)  (with  m  •  Mo)  shows  that  Eq.  (6.63)  holds  quite 
adequately. 

It  is  now  appropriate  to  gather  the  complete  set  of  properly  reduced 
equations  which  determine  Na,  Nb,P,  and  T: 

From  Eq.  (6.16) 


9NA 


-K 


(6.64) 
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From  Eq.  (654) 


^p-5*V2(NBT)  =  K  (6.65a) 

Here,  the  pressure  in  the  form  (from  Eq.  (6.45b)) 

P'ads^T.a^HoiNBT  (6.65a) 

has  been  used.  Most  directly, 

P  =  ali0^NBT  (6.65b) 

Also,  a  diffusion  coefficient  Dn  has  been  defined: 

=  «<*> 

From  Eq.  (6.49) 

Nb9  — ^V(NBT)  (6.67) 

From  Eq.  (658) 

(6.68a) 

where 

DT-jJ^  (6.68b) 

The  development  in  this  section  have  not  depended  on  a  specific  form  for  the 
reaction  rate.  Most  generally. 

K  *  K  (Na,  Nb;  T)  (6.69) 

An  alternate  version  to  Eq.  (6.65a)  is 

2^!1-^DnV2(NbT)  =  TK  +  Nb|[  (6.70) 
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By  Eq.  (6.41a) 


NB3r<<^KelEK  (6.70a) 

at  cv  cv 

Thus  the  last  term  on  the  right  of  Eq.  (6.70)  can  be  neglected  during  propagation. 
The  result  is 


DPV2(NBT)  =  TK  (671) 

ot 

This  is  a  diffusion  equation  for  the  pressure  with  source  =  TK  and  diffusivity 

Dp  =  (Me7a)  A  (672) 

Neither  experiments  nor  the  present  macroscopic  analysis  can  give  an  estimate 
for  Dp.  It  will  be  assumed  that  Dp  is  not  greater  than  Dr;  this  implies  that  not 
much  is  going  on  in  the  cold  material  before  the  approach  of  the  initiation  front. 
This  is  consistent  with  the  previous  assumption  that  the  mass  flow  velocity  can 
not  exceed  the  frontal  velodty.  If  the  flow  resistance  is  sufficiently  great.  Dp 
could  be  considerable  less  than  Dr- 

The  last  topic  to  be  covered  in  this  section  is  the  boundary  and  initial 
conditions  that  must  be  satisfied  by  the  closed  set  of  equations  given  by  Eqs. 
(6.64),  (6.65a),  and  (6.68a)  for  N*,  N&  and  T.  These  are  of  the  form 

(673) 

Only  the  partide  beam  case  (in  cylindrical  symmetry)  will  be  analyzed 
completely.  The  following  sketch  for  $  »  Nq  or  T  is  useful. 


The  initiation  front  is  the  (cylindrical)  surface  over  which  the  temperature  is  Tg 
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and  is  designated  by  the  radial  distance  R(t)  which  is  changing  at  the  rate  dR/dt 
Atr  =  R(t),  Nb  and  T  are  continuous,  however  they  are  expected  to  diminish  very 
sharply  in  regions  II  owing  to  a  drop  in  diffusivity.  The  fluxes  are  continuous  but 
a  discontinuity  in  the  diffusivities  implies  a  discontinuity  in  the  gradients  atr  = 
R(t). 


Since  Eq.  (6.73)  is  of  the  first  order  in  the  time  derivative,  the  initial 

distribution  $(r,o)  (over  all  space)  must  be  specified.  The  initial  time  is  taken 
when  the  materials  in  the  beam  region  is  at  the  initiation  temperature  Tr  and  Nb 
=  Nbr.  No  diffusion  has  occurred  yet  and  the  distributions  are 

T<r'°>={v>a}  *»■> 

NB(r.o)={1^rr<>aa}  (6.74b) 

Also 

Na  (r»  °)  *1  -  Ng{r,  o)  (6.74c) 

Since  Eq.  (6.73)  is  of  the  second  order  in  the  spatial  derivatives,  two  spatial 
boundary  conditions  need  to  be  specified  for  each  region,  conservation  at  the 
boundaries  requires,  for  both  Nb  and  T^ 


ji(o,t)*0  (no  source  at  r*>0)  (6,75a) 

(675b) 

lu  0*0  (no  source  at  infinity)  (6.75c) 
For  the  temperature  at  R(t),  local  thermal  equilibrium  requires 

Ti  (R,  t) »  Tji  (R,  t)  *  Tr  (6.76a) 


This  is  essentially  the  defining  equation  for  R(t),  since  Tr  is  given  as  the  initiation 
temperature.  Similarly  for  Nr, 


NBj(R,t)*NBii(R,t}®NsR  (676b) 

This  last  equality  is  an  assumption  of  constancy  which  is  reasonable  for  the 
quasi-steady  state  behavior  which  is  observed  Also 

NA(fUH-NB(M 

Eqs.  (675)  and  (676)  are  appropriate  sufficient  conditions  to  produce  a  unique 
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solution  for  Na  (r,t),  Nb  (r,t)  and  T  (r,t). 

A  useful  relation  at  the  front  can  be  obtained  as  follows:  Perform  a  volume 
integral  on  Eq.  (6.73)  between  the  limits  R(t)  -  b  and  R  (t)  +  b'  (b,  b'  fixed);  apply 
Gauss'  theorem;  taking  account  that  the  limits  of  integration  are  time  dependent; 
finally  let  b  ->  0  and  b'  The  result  is: 

££  ♦  (r,  t)  rdr  4^1 « (R,  t)  R(t)-  R(t)  j4  (R,  t)=  £  S«(r,  t)  rdr  (6.77) 

It  will  now  be  assumed  that  penetration  of  both  heat  and  gas  into  region  II  is 
small,  thus  the  integrals  in  Eq.  (6.77)  will  be  neglected.  Specifically,  then, 

R|N,(R.t)=-^R|:(NBT)/t„R=NaItR^  (6.78a) 

C„Rj1(R,t)=-CvDTR^/r.li  =  ERRdB  (6.78b) 

where  Er  *Cv  (Tr  -  T0)  is  the  energy  required  to  produce  a  change  in  temperature 
Tr  -  Tq.  Er  is  obtained  from  adiabatic  experiments.  From  the  basic  definition  of 
]n,  as  given  by  Eq.  (6,14),  it  follows 

]Mj(Rt)*NBRVR»NBR^  (&79a) 

and  if  Nqr*  0, 

(6J9b) 

This  completes  the  formal  development  of  the  model  equations.  It  must  be 
emphasized  that  the  basic  material  “point"  is  actually  a  small  macroscopic 
volume  inside  of  which  no  microscopic  modeling  has  been  given.  Caution  must 
be  exercised  in  interpretation.  Since  temperature  measurement,  by  means  of 
thermocouples,  is  the  primary  means  of  obtaining  information,  a  discussion  is 
appropriate. 

It  has  been  assumed  that  the  presence  of  gas  products  is  necessary  to 
account  for  the  observed  rapid  flow  of  heat  into  the  cold  region  of  the  HE. 
Microscopic  heat  transfer  from  the  gas  in  the  porous  channels  to  the  grains  of  HE 
is  likely  governed  by  the  true  thermal  diffusivity  Of  the  HE  grains  for  which  D*  - 
1(H  ctnVsec 

In  Section  3,  the  time  scale  associated  with  initiation  propagation  was  estimated 
be  tj  -  IS  x  1  O'3 sec  in  PBX-9404.  An  associated  thermal  penetration  distance  into 
die  solid  is 
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d(  r  v^Ds  tj  - 10-3  cm 


(6.80) 


Thus,  the  actual  amount  of  HE  material  (not  in  the  beam  region)  that  will  be 
brought  up  to  T  £Tr  will  depend  on  the  morphology  of  the  material,  i.e.,  the 
grain  size  distribution  and  packing  density;  grains  with  d  >  di  will  be 
incompletely  heated.  This  has  bearing  on  the  microscopic  mass  density  of 

reacting  material  which  should  not  be  identified  with  p0/  the  true  mass  density  of 
the  HE  material.  According  to  Eq.  (6.68a),  for  the  quasi-steady  state.  Mo  will  only 

enter  results  through  the  thermal  diffusivity  Dx=  K/p£v;  there  is  thus  no  need  to 
specify  Dr  is  obtained  from  the  experimental  data. 

The  thermocouples  are  in  contact  with  the  hot  gases  and  the  surface  of  the 
solid  grains;  the  measured  temperature  is  therefore  the  temperature  of  the 
porous  channels  through  which  it  is  assumed  that  the  heat  is  being  conducted. 

In  summary  then,  it  has  been  argued,  by  consistent  application  of 
experimental  observations  in  reducing  the  fundamental  hydrodynamic-reaction 
equations,  that  transport  of  mass,  momentum,  and  energy  is  essentially  diffusive 
during  the  observed  initiation  propagation  phase.  The  basic  equations  for  the 


temperature  T  and  the  true  gas  pressure  P  are 

^^-DrVJT.gK 

r  <R 

(6.81a) 

r«R 

(681b) 

and 

i-*v2pa£p&gK 

r<R 

(6.82a) 

-DpR^r/ 

r  =  R 

(682b) 

where 

&«l  +  ||j(NA  +  NB-l)-l 

(683) 

and 

K«K(Na,  N&jD-KtT) 

(684) 
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The  indicated  approximations,  when  appropriate,  render  Eqs.  (6.81)  and  (6.82) 
solvable  without  introducing  Na  and  Nb-  In  general,  the  following  relations, 
involving  Na  and  Nb,  need  to  be  introduced: 

P  =  ^(M?h)N°T  (685) 

~^  =  -K(Na,Nb;T}  (6.86) 

In  the  next  section  Eqs.  (6.81)  and  (6.82)  are  solved  under  appropriately  restricted 
conditions. 
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7.  Solution  of  the  Model  Equation 


In  the  previous  section  equations  were  developed  for  the  determination  of 
the  temperature  T,  the  fraction  of  unreacted  HE  Na,  and  the  fraction  of  reacted 
HE  N&  which  was  allowed  to  be  mobile.  Although  the  autocatalytic  reaction 
model  was  used  as  a  guide,  the  specific  form  of  the  reaction  rate  K(Na,  N&  D 
was  not  specified.  In  Appendix  B  it  is  arguet  lhat  the  autocatalytic  reaction  rate 
is  well  modeled  by 


QKsqRexp 


(7.1) 


in  the  temperature  range  (A  T*  40°  C)  needed  to  analyze  thermal  initiation 
propagation.  In  principle,  if  K(Na,  Nb;  T)  were  known,  the  values  of  qR  and  T« 
could  be  determined;  in  practice  the  values  must  be  determined  from 
experiment.  The  heat  transport  equation,  Eq.  (6.31a)  (with  p  *  p<j)  then  becomes 

2  -  Dr  exp  [{T- Tr]/TJ  (7.2) 

where 

T»  (72a) 

Contributions  from  the  beam  energy  deposition  are  negligible.  Since  no  reaction 
products  are  measured,  Eq.  (7.2)  is  usefully  appropriate.  The  boundary  condition 
at  the  front  is  Eq.  (6.78b): 


■  2vDt  r|~/ r  =  R  =  ,A 


(73) 


It  is  appropriate  to  introduce  dimensionless  variables  into  Eq.  (7.2): 


where 


(7.4a) 


t  -  Evfrr  _  CyT» 


(7.4b) 
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ts  is  the  aduabatic  time  to  explosion.  Also 


when! 


X- 


1 

rs 


rs = 


is  a  diffusion  length.  Additionally 


w  = 


T-Tr 

Ts 


(7.5a) 


(7.5b) 


(7.6) 


Eq.  (7.2)  then  takes  the  form 

(7.7) 


The  boundary  condition  Eq.  (73)  becomes 


-X 


dw 

w 


where 


and 


m 


m 

*TT 


ll-fa 

qR 


(73 a) 


(7.8b) 


(7.8c) 


In  terms  of  the  dimensionless  variables,  Eqs.  (7.7)  and  (7.8a)  have  only  one 
free  parameter: 


t|  Er 

t;  cTi;  (73d) 


Mote  that  this  quantity  doesn't  depend  on  either  qR  or  Dr. 


For  PBX-9404,  adiabatic  experiments  yield 


Er  as  64  cal/ gm  (from  25°  C)  (7.9a) 
Tr  =  553°  K 


Although  values  for  T0  and  Cv  are  obtainable  from  thermal  initiation  propagation 
experiments,  there  is  insufficient  data  to  do  so  in  the  experiments  under  analysis. 
Nominal  values  [Ref.  6],  correct  in  order  of  magnitude,  will  be  used  for 
estimation  purpose.  Thus 


To^265x103*K 
Cv  *  0.4  cal/gm 


These  values  give 


Ts-113e 


(7.9c) 


Similar  estimates  of  ti/t»  for  HBX-1  and  TATB  give  values  of  16  and  23, 
respectively. 


The  quasi-steady  state  equations  are  obtained  by  neglecting  the  time 
derivative  in  Eq.  (77),  This  gives: 


(7.10a) 

_X^/ 

/  X-X  *  dir 

(7.10b) 

in  which 

>1  *1 

(7.10c) 

is  now  the  appropriate  dimensionless  time.  Note  that  Eqs.  (7.10a,  b)  have  no  free 
parameters.  The  quasi-steady  state  solutions  will  be  adequate  as  long  as 


(7.11) 


This  condition  will  be  verified  after  the  solution  to  Eqs.  (7.10a,  b)  is  obtained. 
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In  cylindrical  coordinates  (i.e.,  the  particle  beam  case)  Eq.  (7.10a)  take  the 

form 


(7.12) 


Radial  symmetry  has  been  assumed.  The  general  solution  to  Eq.  (7.12),  obtained 
by  standard  techniques,  is  most  conveniently  expressed  as  follows: 


X(v)  =  e*fr) « 


(7.13a) 


The  quantities  3  and  are  constants  of  integration  which  are  determined  by  the 

spatial  boundary  conditions.  The  gradient  of  w(x)  is 

3W-2(P-1)  4  WXm)2* 

3T  ^  rffe/XmP+i]2  <7*13b> 


This  must  vanish  at  x®  0;  Thus  P  » 1.  In  this  case 


(7.14a) 


8w 

where 

Xq  »  X(x  «  0)  b  (7.14c) 

Xm” 

The  remaining  constant  of  integration  is  obtained  from  the  boundary  condition  at 
the  front;  i.e.,  at  x  «  X  where  w(X) «  0,  thus 


X(X}«1 


(7.15a) 


Since  the  frontal  position  X  is  varying  in  time,  so  is  in  this  manner  the  solution 
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w  becomes  time  dependent  is  a  very  specific  way  which  is  the  essence  of  the 
quasi-steady  state. 


The  solution  for  X*,  as  a  function  of  X  can  be  expressed  in  the  following 
convenient  forms. 


or 


1 

tx2+1 


1±>/i~X2/2 

2 


Xo(X)=4 


XV2  j 


(7.16a) 


(7.16b) 


The  fact  that  Xo  must  be  real  places  a  crucial  restriction  on  X2/^  namely. 


(7.17a) 


With  the  upper  sign  in  Eq.  (7.16b),  the  value  of  X<>  increases  when  X  increases, 

whereas,  with  the  lower  sign  Xo  decreases  when  X  increases,  dearly  the  upper 
sign  corresponds  to  die  physically  observed  quasi-steady  state  behavior;  thus 


X£4 


(7.17b) 


Using  Eqs.  (7.10b)  and  (7.14b),  the  boundary  condition  for  the  flux 
becomes 


Using  Eqs.  (7.16a,  b),  this  becomes 


(7.18a) 


(7.18b) 


This  is  the  equation  that  governs  the  frontal  motion  as  a  function  of  time.  The 
solution  for  6  is 
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e-e'=[v'i-x2/2+iii(i-vM- x2/2)E 


(7.19a) 


In  terms  of  the  original  variables  this  becomes 


t  - 1'  =  titv'l-RViri  +  ln(l  -  Vl-R2/2r?)t 

(7.19b) 

Eq.  (7.19b)  gives  the  time  for  the  initiation  front  to  traverse  the  distance  between 
R'  and  R.  a  comparison  with  experimentally  observed  traversal  times  can  test  the 
model  and  also  produce  values  of  rs  and  tj.  This  will  be  taken  up  in  the  next 
section. 

The  limiting  value  of  R,  from  Eq.  (7.18a)  is 


K  rs  (720) 

Since  it  is  expected  that  initiation  propagation  will  begin  in  the  vicinity  of  the 
beam  edge  (i.e.,  the  collimation  radius),  it  follows  that 

a<R*  (7,21) 

is  a  necessary  condition  for  propagation.  If  the  beam  radius  is  greater  than  R*  an 
explosion  will  develop  on  time  scale  t,-l  msec. 

The  simple  model  presented  in  Section  5  is  obtained  by  letting  Dr  (and 
thus  rj  increase  indefinitely  such  that 


(722) 


Eq.  (7.19b)  then  becomes,  to  first  order  in  1  /r#2, 


R2~K>2 

r»2 


(723) 


The  first  term  in  the  expansion,  with  R'  »  a,  t‘ ‘  to,  agrees  with  Eq.  (5.6)  of  the 
simple  model. 

Attention  will  now  turn  to  the  conditions  for  which  the  quasi-steady  state 
approximation  holds.  Eq.  (7.11)  is  conveniently  put  into  the  form 
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t»8  1  cc. 

!T5eX  1 


(724) 


The  time  dependence  of  the  quasi-steady  state  solution  comes  only  through  its 
dependence  on  the  frontal  variable  X.  From  Eq.  (7.14a) 


Eq.  (7.16b)  can  be  rearranged  (with  the  upper  sign )  to  give 

^-4(1  W1-XV2)2 


(725) 


(726a) 


These  equations  and  Eq.  (7.12b)  are  used  to  evaluate  Eq.  (724)  with  the  results 


1  x% 

n  "7rrr"-.~IT  <  <  1 


2  A4/  J2  Vl-X% 


(726b) 


This  must  hold  for  any  fixed  point  x  <  X.  The  condition  is  least  favorable  for  x  * 
0;  thus  only  this  point  needs  to  be  considered.  From  Eq.  (7.9c)  t,/tj«  1/14  for 
PBX-9404;  thus 


c(x/<E)«^ 


(727) 


A  table  for  C  (X/V2)  follows: 


X/</2 

C(X/V2) 

2 

.025 

A 

.038 

.9 

.066 

.95 

.103 

.99 

248  1 

.999 

.797 

The  quasi-steady  state  approximation  is  clearly  justified  over  essentially  all  of  the 
propagation  region.  As  the  initiation  front  approaches  dose  to  the  limiting  radius 
R\  the  temperature  rate  at  the  origin  undergoes  a  very  rapid  increase.  At  this 
point  in  time  the  dominant  time  scale  changes  rapidly  from  t*  to  t# « ti  and  an 
explosion  develops.  This  warrants  referring  to  R*  as  the  explosion  radius.  These 
conclusions  depend  on  an  assumed  value  for  To;  however,  halving  this  value 
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(and  thus  doubling  ts/ti)  doesn't  greatly  alter  the  conclusions. 
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8.  Model  Parameters  from  Experiments 


The  model  presented  in  the  previous  sections,  in  particular  the  quasi¬ 
steady  state  solution  to  the  heat  transport  equation,  allows  a  comparison  with 
experiment  that  can  produce  the  three  scale  parameters 


(8.1a) 

rs  =  </DTt$ 

(8.1  b) 

jr* 

ii 

(8.1  c) 

The  time  scale  parameter  for  the  non-steady  state. 


ts- 


(8  2) 


can  then  be  obtained  directly. 

Comparison  of  the  observed  frontal  motion  with 


t-t  »t| 


(83) 


for  at  least  three  positions  (two  difference)  will  yield  values  of  t|  and  r9.  If  the 
explosion  radius  R*  is  ascertained,  then 


T*m$[  (8.4) 

is  determined  directly.  Once  ti  and  r,  are  determined,  Eqs.  (8.3)  gives  R  as  an 
implicit  function  of  time.  According  to  Eqs.  (7.19  a,  b),  Xq  and  X  are  then  known 
as  a  function  of  time  and  position.  It  then  follows  that 

IlLjj.T.lR  a  w  » In  X  (83) 

can  be  calculated;  comparison  with  temperature  data  will  then  yield  a  value  for 
T,. 


46 


A  proper  test  of  the  model/  and  the  extraction  of  accurate  model 
parameters  requires  sufficient  amounts  of  data  to  make  any  conclusion 
statistically  significant  As  was  mentioned  previously,  the  available  data  on 
thermal  initiation  propagation  is  taken  from  exploratory  experiments  performed 
prior  to  tile  present  analysis.  The  most  quantitative  detail  is  obtained  from  one 
laser  beam  run  and  two  particle  beams  runs,  all  on  PBX-9404. 

The  laser  case  has  the  best  temporal  resolution  and  shows  details  of  the 
frontal  structure  (see  fig.  A.  2,  lower  left).  In  particular  the  frontal  temperature, 
indicated  by  the  break  in  the  temperature  slope,  agrees  with  the  value  Tr  =  230° 

C  obtained  in  the  accurate  adiabatic  experiments.  The  temperature  was 
monitored  only  at  one  position  in  addition  to  the  metal  -  HE  interface.  Time  is 
measured  from  when  the  laser  was  turned  on.  Although  this  information  is 
conceptually  very  valuable,  it  is  insufficient  for  extracting  model  parameters.  For 
this  reason  the  solution  to  the  model  equations  for  planar  geometry  although 
obtained,  are  not  presented  in  this  report 

The  particle  beam  data  (e.g.  Fig.  A.l)  has  poor  temporal  resolution; 
however  sufficient  information  is  available  to  estimate  ti  and  rs.  Time  of  arrival  of 
the  front  at  two  spatial  positions  can  be  estimated.  The  time  is  relative  to  when 
initiation  occurs  at  the  center  of  the  beam.  The  pertinent  information  for  the  two 
runs  is  abstracted  is  the  following  table: 

acou  8  0.238  (collimation  radius) 


Position 

t(nisec)  f 

I 

a 

.476 

22.2 

19.4 

20  1 

.952 

38.9 

(explosion) 

38*  1 

Table  8.1 

The  accuracy  of  the  time  values  is  judged  to  be  about  30%.  In  run  11  an  explosion 
occurred  before  the  thermocouple  at  position  2  indicated  thermal  initiation.  Since 
initiation  was  indicated  at  this  position  in  run  I  it  is  judged  that  the  average 
explosion  radius  occurs  in  the  vidnity  of  R&  i-e* 

R*  -  .95  cm  (8.6) 


this  immediately  gives 


r* 66  ^  J<i%  m  *673  cm 


(8.7) 
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Defining 


i-£+ta 


and  using  Eq.  (83), 


and  since 


it  follows  that 


tl-  fr-ti 

'  FW-FpC,) 

F(x*a)=o, 

F(Xl)=F(it)°~U44 

tj  •  ti » 18  msec  (nominal) 

tj*  15.7  msec 


(8.8) 

(8.9) 

(8.10a) 

(8.10b) 

(8.10c) 

(an) 


Now  that  r,  and  h  have  been  determined,  Eq.  (83)  can  be  used  to  obtain  the  time 
difference  of  the  arrival  of  the  front  at  two  positions  in  the  propagation  region. 


Quasi-steady  state  propagation  originates  in  the  vicinity  of  the  beam  edge 
at  a  position  which  will  be  denoted  by  a^fj  furthermore  quasi-steady  state 
propagation  is  established,  presumably,  on  a  time  scale  U  « 1  msec  after  initiation 
has  begun.  It  follows  that  foe  absolute  times  in  table  8.1  are  essentially  the  time 
differences  from  to  foe  indicated  position.  These  considerations  allow  a 
determination  of  a^.  Using  foe  information  at  position  2  for  which  HX 2) 113  0,  Eq. 
(83)  gives 


Thus 


a&f 

a? 

=245  cm 


*.133 


The  information  at  position  1,  of  necessity,  gives  the  same  result  In  vie  w  of  foe 
uncertainties  in  foe  data,  the  small  (4%)  difference  in  a<&  and  a<ou  is  perhaps 
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fortuitous;  however  it  does  indicate  consistency  in  the  model  application. 

The  value  of  Ts  cannot  be  obtained  from  the  particle  beam  data  since  the 
latter  does  not  give  the  temporal  development  of  the  temperature  above  Tr.  The 
laser  beam  data  gives  this  information,  but  the  frontal  motion  is  not  known  with 
much  certainty.  For  further  application  the  values  of  Ts  or  tg/ti  given  by  Eq.  (7.9c) 
will  be  used.  An  uncertainty  by  a  factor  of  two  is  judged  to  be  conservative. 

Values  of  qR  and  Dr  can  now  be  estimate  as  follows: 

£ 

4r = if = 4  x  1q3  cal/gm-sec  (8.12) 

Dy  =  ~-  ~  =  0.4  x  103  cm2/sec  (8.13) 

*5  M  lS 

The  magnitude  of  qR  justifies  neglect  of  the  beam  term  (q*b=  10  cal/gm-sec);  it 
also  indicates  that  radical  change  in  the  reaction  has  occurred  at  Tr,  the  initiation 
temperature.  The  large  magnitude  of  the  thermal  diffusivity  Dr  indicated  that 
normal  condensed  matter  thermal  conductivity  is  not  taking  place  (For 
Aluminum  Dr  *  1  cxnVsec).  It's  explanation  is  beyond  the  macroscopic  analysis 
present  in  this  study. 
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9,  Conclusions 


Many  features  of  the  thermal  initiation  propagation  experiments 
described  in  Appendix  A.1  and  Reference  2  can  be  accounted  for  through  the 
diffusion-reaction  equation  (7.2)  with  the  dynamical  boundary  conditions  given 
by  Eq.  (7.3) .  The  quasi-steady  state  assumption  gives  rise  to  the  equation  (7.20): 


t  - 1  =  t^l-(R/R')2  +  ln(l  -  V'l-lR/R'rt  (gl) 

which  describes  the  propagation  of  the  thermal  initiation  front  quite  adequately 
until  the  position  R*  is  readied. 

Detailed  quantitative  validation  of  the  front  propagation  equation  has  not 
been  given  owing  to  the  sparseness  of  useful  data.  There  is  suffident  qualitative 
and  order  of  magnitude  experimental  information  to  allow  a  reduction  of  the 
general  hydrodynamical-reaction  equations  to  the  basic  equations  (7.2)  and  (73). 

The  existing  quantitative  data  for  PBX-9404  allows  for  an  approximate 
determination  of  tj  and  R*;  these  in  turn  yield  a  value  for  the  exothermic  rate  and, 
with  an  order  of  magnitude  assumption  for  the  effective  activation  temperature, 
the  effective  diffusivity,  Dr.  If  reliable  temperature  readings  behind  the  front 
were  available,  a  value  of  the  effective  activation  temperature,  To,  could  be 
obtained  directly. 

The  values 


qg  ■  4  x  103  cal/gm-sec 
Dr  w  0.4  x  103  cm-/sec 


(93) 


do  not  resemble  the  values  inferred  from  Ref.  (6)  which  are 


qK  ■»  39  cai/gm-sec 
Df  *9.7  x  10-4  cmVsec 


(93) 


The  values  are  obtained  as  follows  (with  Tr=  5S3°K) 


qg  *  Z  Q  expC-Ta/Tg} 
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Dt  *  K/p  Cv 


where 


Z  *=  5  x  1019  sec1,  Q  «  500  cal/gm,  Ta  =  265  x  lO^K 
p  =s  1.81  gm/cm3,  cv = 0.4  cal/gm-deg 
k  =  7x  10*4  cal/cm-sec-deg 

These  data  were  obtained  from  experiments  (DSC,  etc)  performed  below  the 
initiation  temperature  Tr.  A 10%  reduction  in  Ta  would  bring  the  values  ofqR  in 
(93)  into  agreement  with  the  values  in  (9.2).  The  large  value  of  Dt  in  (9.2) 
indicates  a  wholly  different  mechanism  for  heat  transport  above  Tr. 

The  principal  applications  of  thermal  initiation  theory  lie  in  the  areas  of 
Directed  Energy  weapons  (particle  beams  and  high  energy  lasers)  and  in  high 
explosives  hazard  analysis.  The  propagation  phenomena  reported  on  and 
analyzed  in  this  paper  are  important  when  the  dimensions  of  the  region  which 
first  experiences  thermal  initiation  is  less  than  the  explosion  radius  R*  The  theory 
derived  here  applies  once  propagation  starts.  Not  enough  information  is 
available  to  determined  all  the  conditions  that  must  be  met  before  propagation 
begins.  A  minimum  stee  region  having  reached  the  thermal  initiation 
temperature  seems  to  be  necessary. 

Further  experiments,  guided  by  die  present  analysis  and  the  work 
described  in  Appendix  A.l  and  Ref.(2),  would  greatly  enhance  the  understanding 
of  thermal  initiation  propagation. 
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Appendix  A.l.  Early  NRL  Propagation  Experiment” 

In  1980,  a  series  of  experiments  were  performed  at  v  $RL  Linac  in  which 

confined  disks  of  HE,  1  inch  diameter  x  1  /4  inch  thick  (about  6  gm),  wc»  s 
exposed  to  non-uniform  irradiation  by  a  32  -  MeV  electron  beam  with  the 
following  pulse  parameters:  0.3  Amp,  1  psec,  and  360  pulses/sec.  this  produced  a 
heating  rate  of  about  9  cal/gm-sec,  and  resulted  in  thermally-induoc  3  e*  L  iosions 
in  6  to  15  sec.  The  HE  samples  were  placed  in  A1  confinement  cells,  sealed  with 
O-rings.  The  beam  entrance  window  was  3/32  inch  thick,  and  the  exit  window 
was  only  1/32  inch  thick,  so  that  blast  fragments  should  move  away  from  the 
accelerator.  The  beam  was  collimated  to  a  diam&er  of  3/16  inch,  so  it  irradiates 
and  heats  only  the  central  region  of  the  HE.  Three  iron-constantan  thermocouples 
(0.001-inch  diameter)  are  imbedded  in  the  rear  face  of  the  HE:  one  at  the  center 
(in  the  beam  area),  and  two  outside  the  beam  arf%  at  radial  distances  of  3/16 
inch  and  3/8  inch  from  the  center.  The  thermoctiipie  response  time  is  about  2 
msec  Two  tests  were  performed  on  each  of  thefollowing  materials:  TATB,  HBX- 
Kcontains  TNT,  RDX  and  Al),  and  PBX-9404  (Contains  94%  HMX).  Beam 
irradiation  was  continuous  until  explosion  occurred. 

In  the  TATB  samples,  explosions  were  quite  weak,  with  no  reaction 
spreading  beyond  the  beam  area  The  HBX-1  samples  exploded  with  more 
violence,  consuming  a  large  portion  of  HE  beyond  the  beam  area.  The 
thermocouple  data  indicate  a  slow  reaction  propagation  velocity,  of  the  order  o*f 
1  cm/sec.  The  PBX-9404  samples  exploded  very  violently,  ripping  up  the  3/4 
inch  thick  aluminum  back  plate  of  die  confinement  cell  into  chunks  of  twisted 
metal.  AU  of  the  HE  was  consumed  rapidly.  Temperature  vs.  time  data  (from  the 
first  run)  from  the  three  thermocouples  near  Initiation  are  shown  in  Fig.  A-l. 
Since  the  thermocouples  were  only  a  millimeter  from  the  highly  conducting  rear 
face  of  the  holder,  they  do  not  indicate  the  tore  temperature  of  the  interior  of  the 
HE.  They  do  indicate  correctly  the  time  of  arrival  of  an  initiation  front  indic&SBd 
by  an  extremely  rapid  temperature  rise.  Anafij of  the  two  PBX-9404  runs  yield 
the  following  data: 


acoii  =  0.238  (collimation  radius) 


Position 

R(cm) 

t(msec) 

I 

E 

Nominal  | 

1 

.476 

22.2 

194 

20 

2 

.952 

38.9 

(explosion) 

38* 

Table  A.1 
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The  accuracy  of  the  time  values  is  judged  to  be  about  30%.  In  run  II  an  explosion 
occurred  before  the  thermocouple  at  position  2  (the  edge)  indicated  thermal 
initjatioti.  Since  initiation  was  indicated  at  this  position  in  run  I,  it  is  assumed 
that  the  explosion  radius  (as  defined  in  the  main  text)  occurs  in  the  neighborhood 
of  position  2.  The  data  indicated  a  propagation  speed  of  about  25-35  cm/sec. 

Recent  data  (1990,  Ref  (2)),  also  obtained  at  the  NRL  Linac,  do  not  indicate 
a  time  difference  between  initiation  at  the  center  and  the  outer  regions.  However, 
these  data  were  obtained  under  rather  different  conditions:  a  wider  beam  (1  /2 
ir  "h  diameter),  a  slower  be*m  heating  rate  (3  cal/gm-sec),  and  using 
thermocouples  with  a  longer  response  time  (40  msecs).  These  experiments  were 
p  tioitned  without  the  benefits  of  the  analysis  presented  in  this  paper.  More 
p<?*tineni  remits  would  require  a  narrower  beam  (£1/8  inch  diameter),  a  shorter 
tinne  response  in  the  thermocouple  ( £  2  msec).  Better  time  resolution  in  time 
recording  is  also  indicated.  Finally,  the  thermocouple  should  be  buried  deeper 
Ubout  half  way)  into  the  interior  of  the  HE  so  that  a  true  profile  of  the 
^  vmpqgatiQn  temperature  can  be  obtained. 


In  hfo  member,  1983,  thermal  initiation  experiments  were  performed  with 
hw  Repetively  Pulsed  Chemical  Laser  (RPCL)  at  the  TRW  facility  in  Capistrano, 
CA  Thr  experimental  arrangement  was  as  follows: 


TT»ehmocoupl«s 

/  "j  * 


Steal 


The  beam  heats  the  outer  face  of  the  steel  plate.  Thermal  conduction 
^  tit  steel  "plate  then  produ  ces  a  temperature  rise  at  the  interface  between 
Lau  FOL  O'vie  thermrcouple  (the  same  type  as  used  in  the  Linac 
vea*  "placed  at  the  kwsrface,  the  other  was  placed  in  the  middle  of 

tea'i^isfv.Te-^nxe  response  is  depicted  in  Fig.  A.2.  the  interface 
'  *%:  >er At'on?  &»#!!&)  rises  -lKK)°c/sec  till  the  HE  melts  in  the  vicinity  of 
-v/flfcc  ,  *  temperature  rise  then  occurs  up  to  560°C,  at  which 

t  '^int  *  v  sharp  temperature1  rise,  indicative  of  intense  exothermic  activity, 
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Figure 
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The  temperature  in  the  middle  of  the  HE  ( lower  left  in  Fig.  A-2)  remains 
unchanged  (~21°C)  until  after  the  above  noted  exothermic  appears.  It  then 
undergoes  a  very  rapid  (<  10  msec)  rise  to  the  vicinity  of  the  thermal  initiation 
threshold  (280°  ±  5°C).  At  this  point  the  temperature  rate  is  considerably 
reduced.  The  precise  structure  of  the  temperature  profile  would  require 
accounting  for  the  finite  (-2  msec)  response  time  of  the  thermocouple.  However, 
the  profile  is  consistent  with  the  propagating  front  structure  described  in  the 
main  text. 

Quantitative  analysis  of  the  frontal  motion  cannot  be  made  since  it  is  not 
known  when  and  where  a  steady  state  initiation  front  began.  This  would 
required  information  from  one  or  more  additional  measurements  in  the  interior 
of  the  HE.  Various  assumptions  result  in  an  estimated  frontal  speed  of  20*90 
cm/sec.  The  order  of  magnitude  is  consistent  with  the  particle  beam  result. 

The  laser  experiment  gives  information  on  what  occurs  when  only  one 
surface  of  HE  is  rapidly  heated.  The  observation  that  the  interface  region  of  the 
HE  reaches  a  temperature  considerably  above  the  thermal  initiation  threshold 
(Tr  «  280°  ±  5°)  shows  that  temperature  alone  does  not  determine  the 
precipitation  of  initiation  propagation.  It  is  probable  that  a  certain  threshold 
volume  has  to  be  brought  up  to  Tr  or  beyond  before  propagation  commences.  In 
the  particle  beam  experiments,  which  produced  rapid  volumetric  heating,  the 
threshold  volume  was  apparently  exceeded.  Quantitative  information  in  this 
matter  is  needed  to  complete  the  understanding  of  the  propagation  phenomenon. 
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Appendix  A.2.  Reaction  Model  Equivalencies 


It  will  be  assumed  that  the  exothermic  heating  rate  can  be  represented  as  a 
function  only  of  the  temperature: 


Q  K  =  qR  G  CD  (A.2-1) 

in  which  qR  is  the  exothermic  heating  rate  at  the  temperature  Tr.  In  general,  since 
G(Tr)*=1, 


lnG(T)  =  ai  (T-Tr)  +  l/2a2(T-TR)2  + 


where 


a2  = 


d2ln  G 
dT2 


In  the  simplest  case 


K(A)aZAe4 


(A2-2) 

(A2-3) 


(A2-4) 


for  which 


Thus,  in  general,  when  the  first  term  in  (A2 -2)  is  an  adequate  approximation. 


(A.2-6) 


can  be  construed  as  an  effective  activation  temperature. 

A  necessary  condition  for  the  adequacy  of  the  approximation  is 

l/2|a2l(T-TR)2«l 

Consideration  of  the  autocataiytic  reaction  rate 

K  *  N  A  <ZA  e  ZB  NBe -$) 


(A 2-7) 


(A2-8) 
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requires  knowledge  of  the  temperature  dependence  of  Na  and  Nb.  In  adiabatic 
uniform  self-heating. 


dNA  K  dNB_r 
dt  —K/  dt  ~K 


cv 


dT 

dt 


=  QK 


The  first  two  equations  are  equivalent  to 


d  Na  _  d  Nb  cv 
dT  "  dT  =Q 


(A.  2 -9) 


(A.  2 -10) 


which  have  the  solution 


N*  Nrr  Cyff-M 

Nar'1‘"Nar  Er  ) 

Nr  !.  .  C»(T-TrN  (A.2-11) 

N^T* . Ir  J 


where  Er  ■  NbrQ  has  been  used.  For  PBX-9404,  cv  -.4  cal/gm-deg,  Er  »  64 
cal/gm,  then  for  Nqr/Nar  “  1/10  and  T  -  Tr  «  40°  C, 

MtM.,25,  Naag.-W.jB 

Er  Nar  Er 

Exponentiation  is  achieved  with  the  relation 


1 4-x®eln(i+*)ae*~f 


(A2-12) 


For  Na/NAr  the  first  term  in  the  exponent  is  an  adequate  approximation  for 
T  -  Tr  S  40°C;  its  only  effect  is  to  make  a  contribution  to  Ta  in  the  amount 


dln-Jr^* 
Nar  / 

dT  / 


-  .19  x  103*  K 


(A2-13) 


with  Tr  *  5S3°K.  Since  Ta  -  263  x  lO^K,  the  contribution  from  Na  is  quite  small 
Henceforth  the  variation  in  Na  will  be  neglected  in  both  adiabatic  self  heating 
and  in  the  Q5S  propagation  process. 
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In  QSS  propagation  the  defining  equations  for  Nb  and  T  are 


-DnV2{NbT)  =  K 
-DTV2T  =  giC 


which  give 

Dfvj  Nb  T  -  Dy  T  =  const 
or 


NB 

NBR  1  +  ^-(T-Tr) 


(A.2-14) 


(A.2-15) 


Although  no  knowledge  of  Dn  is  at  hand#  the  assumption  that  thermal  diffusion 
is  intimately  dependent  on  gaseous  diffusion  leads  to  the  conclusion  that  Dn  and 
Dr  have  the  same  order  of  magnitude.  Thus  the  variation  of  Nb  with  temperature 
does  not  differ  significantly  in  magnitude  between  adiabatic  self  heating  and  QSS 
propagation.  It  is  adequate  to  assume 

NE.Set(T-Tri4(T-TJ1 

Nbr 

where,  for  purpose  of  estimation#  £  *  • 

The  function  GCD  is  conveniently  put  into  the  form 
G(D  -  k  exp  -  3a) + (1  -  k)  exp  (ijg  -  ^ + In 

where 

k  «  ZA  exp  ^ j  |zA  exp  -  ^  +  2b  Nbr  exp  -  ^  j 

(and  Na/Nar  is  set  to  unity).  The  ratio  of  the  first  to  second  term  of  G(T)  at  T  » 
Tp  is  given  by 


(AJ2  - 16) 


(A.2-17) 


(A2-18) 
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thus 


and 


(A.2-19) 


y-  _  1- 
1+XR 


OSxSoo,  l^k^o  (A.  2  -  20) 

On  applying  Eq.  (A.2  -  3),  the  following  basic  relations  are  obtained: 

Ta  =  kTA  +  (l-k)(TB  +  5TRJ)  (A.2 -21) 


k(l-k) 

[Tb-Ta 

2 

Ta 

+  $T, 


(A.2 -22) 


where  T$  =  TR2/Ta.  For  PBX-9404,  Tr  =  553%  Ta  «  26.5  x  103%  using  §  *.0063  it 

follows  that  T,  *  11.54°,  £TS » .073,  §Tr/2  ~  1.74.  The  £  dependent  part  of  the  first 
term  in  (A~2  -  22)  is  seen  to  be  negligible.  Thus 


.{-4a +^[3^A+073f}(I^)a 


(A.2  -  23) 


Without  specific  knowledge  of  xR  and  Tb  -  Ta  it  is  necessary  to  consider 
extremal  cases.  At  the  extremes  k  «  0, 1  (xR  «  0), 

EXTREME  gj  =  -  -021  (A.2 -24) 

Kite  coefficient  a2  is  maximized  at  k*  1/2  (xr  =  1) 

MAXIMUM  g2«[~  .021+1^6^  +  0.73^1^)*  (A3 -25) 


The  maximum  exceeds  the  magnitude  of  the  extreme  value  only  when 

2&~£<-.635  {Tb < Ta) 

or  “  (A3 -26) 

Jirl£ > 507  (Tb>Ta) 

•a 
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Thus  for  a  considerable  range  of  values  of  Tb  -  Ta  and  for  all  values  of  k  (or  xr)/ 


I|a2|p-^S)2 <  '021  pT^f  <<1  (A2-27) 

can  be  taken  as  a  rather  conservative  condition  for  die  validity  of  the  first 
approximation. 

In  QSS  propagation  the  maximum  value  of  CT  *  Tr)/Ts  is  In  4  =  1386  (AT  s 
16*Q,  thus 

£ia2lP^k)2<  -04 
More  generally,  when  T  -  Tr  £  40°C, 

<-25 

This  is  probably  the  limit  for  which  the  approximation  is  useful. 

In  summary,  the  analysis  shows  (1)  that  the  simple  form 

QK-qge&tr-T‘> 

can  be  deduced  from  more  complex  forms  and  used  for  a  limited  temperature 
range  and  (2)  the  analysis  of  data  for  a  limited  temperature  range  AT  £  40°C  can 
only  yield  two  reaction  parameters,  qR  and  To)  more  complex  modeling  is 
inappropriate. 
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